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ABSTRACT 

We present new Planetary Nebula Spectrograph observations of the ordinary ellip- 
tical galaxy NGC 4494, resulting in positions and velocities of 255 PNe out to 7 
effective radii (25 kpc). We also present new wide-field surface photometry from 
MMT/Megacam, and long-slit stellar kinematics from VLT/FORS2. The spatial and 
kinematical distributions of the PNe agree with the field stars in the region of overlap. 
The mean rotation is relatively low, with a possible kinematic axis twist outside 1 R e . 
The velocity dispersion profile declines with radius, though not very steeply, down to 
~ 70 km s _1 at the last data point. 

We have constructed spherical dynamical models of the system, including Jeans 
analyses with multi-component ACDM-motivated galaxies as well as logarithmic po- 
tentials. These models include special attention to orbital anisotropy, which we con- 
strain using fourth-order velocity moments. Given several different sets of modelling 
methods and assumptions, we find consistent results for the mass profile within the 
radial range constrained by the data. Some dark matter (DM) is required by the data; 
our best-fit solution has a radially anisotropic stellar halo, a plausible stellar mass- 
to-light ratio, and a DM halo with an unexpectedly low central density. We find that 
this result does not substantially change with a flattened axisymmetric model. 

Taken together with other results for galaxy halo masses, we find suggestions 
for a puzzling pattern wherein most intermediate-luminosity galaxies have very low 
concentration halos, while some high-mass ellipticals have very high concentrations. 
We discuss some possible implications of these results for DM and galaxy formation. 

Key words: galaxies: elliptical and lenticular — galaxies: kinematics and dynamics 
— galaxies: structure — galaxies: individual: NGC 4494 — dark matter — planetary 
nebulae: general 
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1 INTRODUCTION 



The idea that some sort of dark matter (DM) dominates the 
mass of the universe arose primarily in the 1970s fr om obser- 
vatio n s of the outer kinematics of spiral galaxies dFreemanl 
19701; lOstriker fc Peebles! 1 19731 ; iRubin et all 1 19781 ; iBosmal 



20031; lO'Sullivan fc Ponmanll2004l; IPellegrini fc Ciottil 120061 : 



1981), and is now we ll established on cosmological grounds 



I Hinshaw et al.ll2008h . As part of the "concordance" ACDM 
model — including cold dark matter (CDM) and a cosmolog- 
ical constant (A) — DM has been invaluable in describing the 
formation, evolution, and stability of cosmic structures. 

Yet despite the successes of the concordance model on 
cosmological scales, on the scales of galaxies and galaxy clus- 
ters, several observational discrepancies continue to chal- 
lenge cosmological complacency. For instance, ACDM sim- 
ulations of th e formation o f DM haloe s predict a steep 
central cusp dNayarro et all Il99fj , 1 1991 hereafter NFW; 
M oore et al. 

I Il999h . a t variance with many observations 
of late-type galax i es (Ide Blok. Bosma. fc McGaugbl | 2003 |; 
Gentile et al|2004l ; iGentile et alj|2005l; iGilmore et al.ll2007l; 
Salucci et al.ll2007l ; ISpano et a l. 2008; K uzio de Narav et al.l 
20081 ). It remains to be seen whether this discrepancy can 



be traced to observational problems, to oversimplified pre- 
dictions of DM halo properties (concerning in particular the 
inclusion of baryonic effects), or to a failure of the ACDM 
paradigm. 

The default expectation for baryons is that their dissi- 
pation would exacerbate the central DM density problem, 
but it is possible that processes such as feedback in a com- 
plete galaxy formation picture would decrease the DM den- 
sities. 

Given the long-standing DM puzzles in late-type galax- 
ies, it is increasingly important to examine the mass content 
of early-type galaxies (ETGs), i.e. ellipticals and lenticulars. 
Since these objects are usually free of cold gas and have their 
stars moving in random directions, their kinematics are more 
difficult to sample than in disc galaxies and their dynamics 
harder to model. 

The first obstacle comes from observations. Outside 
~ 1-2 R c (where the projected effective radius R e is that 
encircling half the total light of the galaxy), the decreasing 
surface br ightness makes kinematics measurements very dif- 
ficult re.g.lHallidav et ai]l200ll ; ISamurovic fc Danzigerll2005l ; 



ICorsini et al. 20081 ) . DM constraints from this technique 
are therefore possible for only a fraction o f galaxies (e.g. 



Kronawitter et aljlioool ; iGerhard et~aT]|200ll ; iThomas etail 
20071 ). 



Globular clusters have classically been used as mass 
tracers in the haloes of "bright" galaxie s, but the kine- 
matic al samples in "ordinary" ETGs (|Capaccioli et all 
1992) have so far been too small for strong constraints 
(e.g. iPuzia et al] |2004 iBergond et~aH l200rj) . Similarly, 
mass studies using X-ray emission are biassed toward 
massive, group or cluster-central galaxies, while ordinary 
ETGs are much more difficult to probe (see iPaolillo et al.l 
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Humphrey et al.ll2006l ; IPellegrini et al.ll2007l ). 

Hence it is not surprising that there is no unbiassed, 
systematic survey of the detailed DM properti es of ETGs, 
comp arable to what is available for late-types (I Persic et al.l 
Il999) . 

A way forward is offered by planetary nebulae (PNe), 
which have become established probes of the stellar 
kinematics in ETG haloes, where their bright 5007 A 
[O in] line stands out against t he faint galaxy back- 
ground (e.g. Cia-rdullo et al. I Il993l; lArnaboldi et all Il996l; 
Arnaboldi et al] ll99Sl; iNapolitano et al] |2002|; IPeng et al' 



2004 ). In addition to their use as large-radius mass probes, 



PNe have the additional advantage of tracing the kinemat- 
ics of the underlying field stars, and so their orbital proper- 
ties can be used to constrain the formation history of the 
bulk of the host galaxies' stars. In the last few years, a 
novel observational techniqu e of "counter-dispersed imag- 
ing" (|Douglas fc Taviorl lT999) has overtaken the traditional 
multi-object spectroscopic follow-up of PNe candidates dis- 
covered in narrow-band images. This efficient approach to 
collecting large samples of PN velocities relies on the com- 
parison of pairs of slitless spectroscopic images taken with 
the dispersion element rotated by 180° between the expo- 
sures, and has bee n successfully implemented with several 
different variations jPouglas et al]|200fj| ; IMendez et"alll200ll ; 
iTeodorescu et al.ll2005l) . 

The Planeta ry Nebula Spectrograph (PN.S; 
iDouglas et aill2002l ) is the first custom-designed instrument 
for counter-dispersed imaging. This instrument's two spec- 
trographs cameras, fed by a beam-splitter arrangement 
at the grating module, operate simultaneously, securing 
image pairs from the same collimated beam and via the 
same filter at the same temperature. The PN.S has begun 
producing large k inematical samples o f PNe in a variety 
of ga laxy types dMerrett et al] 120061; iRomanowskv et all 
20031. hereafter R+q3~lDouglas et al.| l2007i. here after D+07; 



Noordermeer et al] 120081 , ICoccato et al] 120081 ). The first 



results from a small sample of ordinary ellipticals yielded a 
surprisingly pseudo-Keplerian decline in the observed ve- 
locity dispersion profile, more suggestive of galaxies with a 
constant mass-to-light ratio (M/L) than of DM-dom inated 
systems dCiardullo et ail 1 19931 ; IMendez et ail l200ll ; R+03; 
ISluis fc Williamsll2006l ; D+07) . 

These PN observations have sparked several main in- 
terpretations. The first is that the data are c ompatible 
with A CDM expectations dMamon fc L okas 2005 , hereafter 
ML0 5; iDekel et al] 120051 . hereafter D+05; lOnorbe et al] 
I2007D . The second is that the central DM concentra- 
tions of many ordina ry ETGs are unexpectedly low 
|Napolitano et ahl 120051 . hereafter N +05), akin t o sim- 
ilar conclusions am ong late-types ((Kassin et al] 120061 ; 
iMcGaugh et al]|2007l ). The third is that a modified form 
of gravity, rather than DM, dominates the outer regions 
of ga l axies dMilgrom fc Sandersjl2003l ; iBrownstein fc Moffatl 
l200fj ; lTiret et al]|2007l )~ 

The issues raised by D+05 in particular (relating to 
anisotropy, viewing angle, and star-PN connections) were 
explored by R+03 and D+07 in the context of NGC 3379— 
the paragon of an ordinary elliptical, with so far the best 
available kinematical data. They concluded that the ACDM- 
based interpretations were not convincing, although the sit- 
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Table 1. NGC 4494: basic data. 



Table 2. Log of NGC 4494 PN.S observations 



Parameter 



Value 



Reference 



R.A. (J2000) 
Decl. (J2000) 

^sys 

(m - M) 

A B 

my 

TUB 

M v 
M B 
(B- 

Re 
TO 

e 

04 
PA 



V) 



12 h 31 m 24 s 
+25° 46' 30" 
1344 km s" 1 
31.0 mag 
0.09 mag 
9.74 ± 0.10 mag 
10.55 ± O.lOmag 
-21.26 ± O.lOmag 
-20.5 ± O.lOmag 
0.81 mag 
48".2±3".0 
150 km s" 1 
0.162 ± 0.001 
0.16 ±0.07 
-0.9° ±0.3° 



NED 1 

NED 

NED 

Tonrv et al. (2001) 2 

Schlegel et al. (19981 

Sec. 

Sec.lXTl 

Sec. [311 

Sec.O 

Goud frooii et al. (19941 

Sec. [XT] 

HyperLeda 3 

App.fA] 

App.E] 

App.S 



NOTES - (l): |http://nedwww .i pac.caltech.edu/ 

(2) : corrected by -0.16 mag (sec Jensen ct al. 2003) 

(3) : |http://leda.umv-lyonl.fr/| { Paturel et alj|2003li . 



uation may be less cle ar under more general dynamical 
modelling assumptions (|De Lorenzi et aL l2008bl . hereafter 
DL+08b). DL+08b explored flattened galaxy models in de- 
tail, and found a fairly wide range of mass profiles permitted 
for NGC 3379. Although they did not directly test ACDM- 
based mass models, they made an a posteriori comparison 
to the D+05 simulations. They found that the data permit- 
ted total mass models comparable to D+05, but that these 
required stronger radial anisotropy than found in the simu- 
lations. 

With various efforts underway to elucidate the implica- 
tions of the data, there also remains the task of producing 
the observational results for a broad sample of galaxies. Here 
we continue toward this goal by presenting new photomet- 
ric, kinematical, and dynamical results from a second galaxy 
with preliminary results in R+03: NGC 4494. 

NGC 4494 is a LINER El galaxy located in the Coma I 
cloud on the periphery of the Virgo cluster, at a distance 
of 15.8 Mpc and with a total luminosity of Mb ~ —20.5, 
a value close to the L* characteristic luminosity. These and 
other global parameters are listed and referenced in Table [1] 
The stellar component of NGC 4494 is somewhat discy, with 
a steep central cusp and signification rotation, exhibiting few 
peculiarities other than central discs of dust (~ 100 pc) and 
kinematically-decoupled s tars (~ 400 pc), wh i ch may imply 
a rec e nt merging episode dBender et al. 19941 ; ICarollo et al.l 
Il997l ; iTran et all l200ll ; lLauer et al.l |2007| ). X-rav observa- 
tions reveal very little hot gas in the halo, implying either a 
recent interacti on which has depleted the gas, or a very low- 
mass DM halo (|Q'Sullivan fc Ponm an 2004; Fuka zawa et al.l 
l200fj ). 

Our new study of NGC 4494 follows the pattern of re- 
duction and basic analysis of NGC 3379 in D+07, but goes 
further in establishing a stronger framework for deriving and 
interpreting the mass profile. These include the critical as- 
pects of exploiting higher-order velocity information in the 
data to break the mass-anisotropy degeneracy, and fitting a 
suite of physically-motivated multi-component mass models. 

The paper is organized as follows. Section [2] presents a 



Date 


Integration 


Seeing 




[hours] 


[FWHM] 


2002 March 7 


1.3 


1.9" 


2002 March 8 


1.0 


1.6" 


2002 March 9 


2.8 


1.6" 


2002 March 13 


0.7 


2.1" 


2002 March 14 


5.3 


1.4" 


2003 February 28 


0.5 


1.6" 


2003 March 1 


0.5 


1.6" 


2003 March 2 


0.5 


1.6" 


2003 March 3 


1.2 


1.8" 


2003 March 4 


0.5 


1.0" 



standard reduction of the PN.S observations of NGC 4494. 
In Section [3] we present the spatial and kinematical proper- 
ties of the PN system, comparing them to newly acquired 
long-slit stellar data. We analyse the system's dynamics in 
Section [4] and put the mass results in context with other 
galaxies in Section [5] Section [6] summarizes. New surface 
photometry for NGC 4494 is tabulated in Appendix [A] and 
some technical aspects of the dynamics are described in Ap- 
pendices [B] and [C] 



2 PN.S OBSERVATIONS OF NGC 4494 

We present a new set of PN.S data on NGC 4494, starting 
with a description of the observations and basic data reduc- 
tion in £|2 . 1 1 and addressing the completeness and possible 
contamination of the sample in 



2.1 Observations and data reduction 

PN.S observations of NGC 4494 were performed during two 
runs (March 2002 and February-March 2003; Tabled under 
a variety of seeing conditions, for a total integration time of 
~ 14 hours (equivalent to 6.2 hours at a constant seeing of 
FWHM = 1", for a given S/N in the final co-added imagsQ). 
Preliminary results from the 2002 data were presented in 
R+03. 

The systemic velocity of NGC 4494 (see Table [TJ was 
matched with our "B" narrow- band (FWHM= 31 A) filter 
tilted by 6° to achieve a shift of 7.16 A of its nominal cen tral 
wavelength A c = 5033.9 A (see also iDouglas et al. I l2002l for 
further details). The instrument position angle was set to 0°, 
such that the grating dispersed in the North-South direction. 

The observations were processed as detailed in Section 
2 of D+cfl. An important step in the procedure is identify- 
ing PNe in the dispersed images by their pointlike emission, 
since they are spatially and spectrally unresolved objects. 
Potential sources of contamination are noise, foreground 



1 In background-limited observations, 



S/N OC SEEING" 



SO EXP TIME OC SEEING 



S/N oc v'exptime and 
-2 



2 Data reduction is via a dedicated pipeline written in the IRAF 
script language (distributed by NOAO, which is operated by 
AURA Inc., under contract with the National Science Founda- 
tion), with some additional routines written in FORTRAN. 
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Table 3. Catalogue of PNe in NGC 4494 



ID 




R.A. 


Decl. 


Wavelength 


^hclio 


PNS-EPN 


- 


(J2000) 


(J2000) 


[A] 


[km s' 1 ] 


NGC4494 


1 


12:31:45.3 


25:47:26.9 


5029.3 


1349 


NGC4494 


2 


12:31:45.4 


25:45:20.5 


5029.4 


1358 


NGC4494 


3 


12:31:44.1 


25:48:50.0 


5029.8 


1377 


NGC4494 


4 


12:31:43.7 


25:49:07.5 


5029.1 


1341 


NGC4494 


5 


12:31:43.5 


25:45:47.4 


5028.7 


1314 


NGC4494 


6 


12:31:42.8 


25:45:46.2 


5027.6 


1250 


NGC4494 


7 


12:31:42.5 


25:47:02.8 


5028.3 


1290 


NGC4494 


8 


12:31:41.6 


25:48:28.9 


5017.9 


670 


NGC4494 


9 


12:31:41.2 


25:46:04.5 


5028.2 


1275 


NGC4494 


10 


12:31:40.0 


25:47:29.4 


5029.2 


1346 


NGC4494 


11 


12:31:36.3 


25:49:59.1 


5026.8 


1199 


NGC4494 


12 


12:31:35.2 


25:47:56.3 


5030.8 


1437 



NOTES - The PN.S data follow the usual definitions, and the ID is chosen to be compatible with proposed IAU standards 
(EPN standing for "Extragalactic Planetary Nebula"). Velocities include the heliocentric correction of 5.7 km s _1 . 




NGC4494 PN.S/R 
■ m » * • 




Figure 1. Dispersed images of NGC 4494 from the PN.S, with 
close-up of the central regions; the "Left" and "Right" camera 
arm images are marked by PN.S/L and PN.S/R. The galaxy cen- 
tre is marked by X . The images have been median subtracted 
and Gaussian smoothed in order to enhance the emission lines. 
The stellar background precludes practically any detection for 
R < 20", and the crowding is evident. Some sample pair associ- 
ations are framed: in the green box two pairs clearly separated 
in the two arms; in the blue box three emissions in the Right 
arm collapse into two in the Left arm; in the red box there is a 
complicated set of 4 or 5 emissions. 



and background continuum objects, and background galax- 
ies with line emission ([O n] at redshift z ~ 0.35 or Ly-a 
at z ~ 3.14). The latter two categories can often be iden- 
tified by their extended shape in the images, while noise 
is rejected by requiring matching detections in both arms. 
The point-like objec ts are extracted using t he automated 
software SExtractor (I Bert in fc Arnoutdli99rjh . with extrac- 
tion parameters optimized by experimenting on simulated 
data sets. The automated selection is complemented by an 



"eyeball" check on every PN candidate by at least three ob- 
servers. This check also picks up some PNe missed by SEx- 
tractor in certain regions of the image (near star-trails and 
the galaxy centre). 

PN.S observations near galaxy centers are generally 
hampered by the increasing dominance of the stellar back- 
ground, and by crowding as a consequence of the central 
concentration of the PNe themselves. These effects lead 
to reduced detection completeness, as well as to confusion 
in assigning PN pairs detected in the left and right arms, 
and were more significant in NGC 4494 than in NGC 3379 
because of a poorer median seeing, higher stellar surface 
brightness, and larger number of PNe in the case of the for- 
mer (see Figure [1]). Fortunately, the central galactic regions 
are not of much concern to us, since the power of the PN.S 
lies in probing the outer parts of galaxies. We must simply 
keep in mind that the PN-based results inside ~50" should 
be considered unreliable. 

After rejecting ~ 30 objects for lack of unanimous ap- 
proval from the independent observers, we obtain a remark- 
able full sample of 267 PN candidates (Table [3jl. These ob- 
jects may still include a small number of background line- 
emitter contaminants and pair mismatches (see Section r2.2[) . 
Their spatial distribution is shown in Fig. [2] 



2.2 Completeness and contamination 

It is important to understand the completeness of our PN 
sample as a function of radius and magnitude (see £)3.1I) . 
Therefore, we followed the procedure used in D+07 for es- 
timating the fraction of sources missed by our PN search 
algorithm. We added to our images a simulated sample of 
400 point-like sources randomly placed, with a luminosity 
function following that expected for the PNe of NGC 4494 
and co vering a range up to 1 m ag fainter than the bright-end 
cutoff l|Ciardullo et al.lll989a) app ropriate to this galaxy, 
m* ~ 26.0 mag ( Jacobv et al.l 19961 ) . where the [O III] mag- 
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NGC 4494 

100 -100 -200 -300 




-100 



- -200 



-300 



100 - 100 

ARA (arcsec) 

Figure 2. Spatial distribution of objects around NGC 4494. The 
PNe are marked as filled green boxes, with red X symbols show- 
ing objects discarded as phase-space "outliers" (see text). Surface 
brightness isophotes from (Megacam (/-band) are shown by con- 
tours, and the blue solid line marks the boundary where the sur- 
face photometry is unaffected by the bright star to the Northeast. 



nitudes are defined as 

tosoo7 = -2.5 log F 50 o7 - 13.74 



|jacobvlll989l ) . We ran the same SExtractor procedures used 
for real data, finding that ~ 90% of the sources were re- 
covered outside R ~ 3R C (~ 150"). The detectability de- 
creased to ~ 75% at R ~ R e , falling sharply inside this 
radius. Inside R ~ 20", even the brightest (m*) PNe were 
not recovered. In the real data, the completeness was higher 
because of the additional visual checks, so we correct our 
completeness fraction upward by using the ratio of auto- 
matically and visually recovered PNe in the real data, as a 
function of radius. This correction is ~15% and ~10% in 
the inner and outer regions. Overall, our 50% completeness 
limit was (m* + 1.2) ~ 27.2 mag. 

We next try to eliminate any residual false PN iden- 
tifications, which could be due to unresolved background 
emission-line galaxies, or to PN pair mismatches in crowded 
regions. As in D+07, we recognize such cases by their ap- 
parent "velocities" which are outlying from the general dis- 
tribution defined by the bona fide PNe. To this end, we 
use th e "friendless" algorithm introduced in iMerrett et al.l 
|2003l ) which flags objects deviating by more than n x a 
from the velocity distribution of their N nearest neighbours. 
With n = 2.5 and N = 15, we find 12 outliers among the 
NGC 4494 PN candidates (see Fig. [2). 

Almost all of the friendless objects are in areas of the 
field that are either noisy or prone to pair confusion; one 
or two remaining outliers are probably background galaxies. 
An alternative approach calculates an azimuthally-averaged 
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Figure 3. Distribution of line-of-sight velocities of PN candidates 
around NGC 4494, as a function of radius, and relative to the 
systemic velocity (1344 km s — 1 ). Red X symbols mark objects 
designated as outliers by the friendless criterion (one object is off 
the plot at 108", 2518 km s" 1 ), and green boxes show the bona 
fide PNe. The dotted line shows the 3 a velocity envelope. 



projected velocity dispersion profile using a moving-window 
technique, and rejects objects outside the 3 a velocity enve- 
lope (see Fig. [3}. This approach finds the same outliers as 
the friendless algorithm, except in the very central regions 
which are not of importance to our study anyway. The fi- 
nal PN catalogue for NGC 4494 therefore has 255 objects; 
this is 3.5 times the 73 PNe reported in R+03, with most 
of the increase coming from a greatly improved method of 
reducing and stacking the raw data. 

We have compared the R+03 preliminary catalogue 
with our new one and found 70 PNe in common between 
the two datasets, with no velocity offset between the two 
samples (for comparison, in D+07 we found an offset of 
~ 30 km s" 1 for NGC 3379). The distribution of veloc- 
ity differences between the R+03 sample, Vr+03, and the 
present sample, Vpnos, AV = Vk+03 — Vpnos is nearly Gaus- 
sian (skewness ~ —0.4 and kurtosis ~ 2.9) with a standard 
deviation of 18 km s _1 — which is smaller than the typical 
instrumental error of the PN.S (20 km s , see D+07). We 
find no systematic trend of velocity difference with spatial 
position. 

We have also cross-checked the Ijacobv et al. (1996, 
hereafter J+96) photometric-only catalogue of 183 PN can- 
didates, all of them in our field of view. We recover only 
121 of their PNe; of the ones we miss, two-thirds are not 
from the photometrically complete part of their catalogue 
(wi5007 < 27.0, R > 1.0'). Photometric catalogues are known 
to be liable to false d etections if the off-b and images are 
not deep enough (e.g. lAguerri et al J 12005). and the frac- 
tion of the photometric PNe that we have confirmed is 
consistent with typical recovery rate of photometric candi- 
dates in spectroscopic follow-up (e.g. ICiardullo et al.lll993l ; 
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lArnaboldi et al.ll 19981 ) — where the losses can also be due to 
different factors like astrometric problems, fiber positioning 
and losses etc. On the other hand, the vast majority of the 
objects detected only by the PN.S were below the complete- 
ness limit of J+96. 

A comparison of spatial positions shows a large average 
difference in RA and DEC (~ — 2" and ~2" respectively) 
between our positions and the ones from J+96. However, 
the scatter is small (0".6 and 0."5, respectively), so there 
is simply a constant offset between the astrometric systems, 
which will not compromise the kinematical and dynamical 
analyses of this paper. 

As a final check we have compared the [OIII] magni- 
tudes m5oo7 between the two datasets, finding a scatter of 
~ 0.4 mag for the whole matched sample, which reduces to 
~ 0.3 mag for the brightest objects (within 0.3 mag of m*). 
Such scatter is consistent with the typical photometric un- 
certainty of the PN.S, of the order of 0.2-0.3 (see also D+07) , 
once the photometric errors quoted by J+96 are subtracted 
off in quadrature. 



3 PN DISTRIBUTION AND KINEMATICS 

We now present the basic properties of the field stars and 
PNe in NGC 4494, including their distributions in space and 
velocity. Since an important assumption in strengthening 
our models is that the PNe are a fair tracer population of 
the field stars, we compare throughout the properties of the 
stars and PNe. We examine the spatial distribution in ^3. II 
the rotation field in H3.2I the velocity dispersion in H3.3I and 
the kurtosis in ^3.41 



3.1 



Surface photometry and PN spatial 
distribution 



The surface photometry of NGC 4494 available in the lit- 
erature in£hodes_^57]42ased observations in the V and / 
bands (lLauer et al.ll2005 | ). and ground -based CCD observa- 
tions in BVI (iGoudfrooii et~aLlll994l . hereafter G+94). As 
neither study extends far enough in radius to cover the re- 
gion probed by the PNe (to ~ 6'), new surface photome- 
try was necessary. To this end we have use d imaging data 
taken on January 27, 2006, with Megacam jMcLeod et al.l 
|2000| ) mounted on the MMT 6.5m telescopfl as part of a 
programme to image the globular cluster (GC) systems of 
nearby ETGs. A mosaic of 36 2kx4k backside-illuminated 
CCDs provide a 24'x24' field of view with a 0.08" native 
pixel scale. Our data consist of 3x 45 sec exposures of 
NGC 4494 using the SDSS g' filter, read out using 2x2 
binning, with a pixel scale of 0.16". The detailed data reduc- 
tion and surface photometry analysis are beyond the scope 
of this paper, but we provide a brief summary here. 

The initial data reduction (including overscan, trim- 
ming, flat-fielding, and bad pixel identification) was per- 
formed using the standard routines in the IRAF MSCRED 
package. World Coordinate System solutions were then 
found for each exposure by comparison with the 2MASS 



catalogue (jjarrett et al.| [2003+ The images were projected 
to the tangent plane using SWarrJE and stacked in IRAF to 
create the final corrected images with a seeing-limited reso- 
lution of FWHM ~ 1.4". Isophotes from the reduced image 
are shown in Fig. [2j where it is evident that contamination 
from a bright star becomes problematic beyond ~ 2' in the 
Northeast direction. Therefore we have excluded this region 
as marked in the Figure from our surface photometry anal- 
ysis. 

To convert the surface photometry into simple photo- 
metric profiles suitable for modelling, we used the IRAF 
tasks IMSURFIT and ELLIPSE to subtract the background 
light and fit elliptical isophotes. The resulting profiles of sur- 
face brightness, ellipticity, and position angle are shown in 
Fig. U mapping from the (/-band to the V-band with an ar- 
bitrary normalisation required to match the literature data 
over the radial range of R = 1.8"-60". In this Figure, and in 
general in this paper, the radius used is from the projected 
intermediate axis R, n , which is related to the major axis R a 
and ellipticity e by R m = R a (l — e) 1//2 . The data sets are 
seen to be generally consistent, except notably inside ~ 4", 
where the Megacam data are affected by seeing. We there- 
fore combine the literature and Megacam data into a single 
l/-band surface brightness profile /j,v(R), listed in Table |A"T1 
of Appendix lA") together with ellip ticity and isophote shape 
parameter cu (Bender et all ll 988) values. 

In order to make use of the stellar luminosity profile in 
the dynamical models, we parametrize the surface brightness 
profile by the Sersic law: 



M (i?) - K0) oc (i?/a s ) 1/m , 



(1) 



where as is a scale length and m describes the "curvature" 
of the profile l|Sersid[l968l ). Fitting the regions outside the 
inner disc (R > 5"), we find as = 0.115", m = 3.30, 
and no = 14.82 mag arcsec -2 . These parameters trans- 
late into an effective radius of _R C = 49.5"; a more de- 
tailed multi-component fit including the central regions 
yields R c — 48.2 " ± 3.0", which i s similar to other lit- 
eratu re estimates dFaber et al.l 1 19891 ; Ide Vaucouleurs et al.l 
Il99ll ). The total extinction-corrected luminosity in the V- 
band is 2.64 x lO lo Ly, , or M v = -21.26; the uncertain- 
ties in the outer surface brightness profile yield a (model- 
dependent) total luminosity uncertainty of ~5-10%. To es- 
timate the same quantities in the B-band, we use the typical 
de-reddened color of (B - V) = 0.81 from G+94 to find 
L B = 2.37 x 10 10 Lb, q and M B = -20.45. These and other 
global parameters for NGC 4494 are listed in Table [T] 

We next compare the spatial density of the PNe with 
the field stars, using the PN number density complete to 
m* + 1 (Section I2.2|l . Given an arbitrary normalisation, 
the PN profile matches the stellar photometry remarkably 
well (Fig. [5} — -as also genera l ly fo und in a larger sample 
of galaxies by ICoccato et alj l|2008l ). Therefore we do not 
see a pattern emerging to support the scenario proposed 
by D+05, wherein the observed PNe trace a young stellar 
sub-population generated in a merger, with a steeper spatial 
density profile than the bulk of the field stars. 



3 The MMT Observatory is a joint facility of the Smithsonian 
Institution and the University of Arizona. 
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manuel Bcrtin 



software 
and it 



developed 
is publicly 



by Em- 
available at 



|http://terapix.iap.fr/rubrique.php?id-rubrique=49/| 
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Figure 4. Surface photometry of NGC 4494, as a function of 
the intermediate axis r adius. Literature data iGoudfrooi i et all 
11994 lLauer et alj|2005h are shown by red circles, and Megacam 
data by black stars (this paper), with the uncertainties shown by 
error bars. Top large panel: composite V-band surface brightness 
profile. The solid line is the best fitting Sersic (1968) formula for 
R > 5". Bottom small panels: Profiles of ellipticity, position angle 
and fourth-order isophote shape coefficient 04, with G+94 data 
also shown. 



The normalisation of the PN profile corresponds to a 
specific density parameter of PNe per unit l/-band galaxy 
luminosity of a v, 1.0 = (11. 5±3.O)(1O 9 Z/0) _1 . This translates 
to a B-band density of a B ,i.o = (10.3 ± 2.6)(10 9 L©) _1 , and 
over the standard magnitude interval to m* + 2.5, as, 2. 5 = 
(42 ± ll )(10 9 I/p))~ 1 , which i s consistent with the previous 
result of lJacobv etail (Il996l ): a B ,2.5 = (31 ± 3)(10 9 L o ) _1 
(assuming our adopted extinction and distance). The spe- 
cific frequency is fairl y high compared to other ETGs of 
simil ar luminosity (cf. ICiardullo et al.l 120051 . iBuzzoni et al.l 
120061 ). The 20% discrepancy between the PNLF and SBF 
distances to NGC 4494 will be investigated in a future pa- 
per on PNLFs. 



3.2 Rotation 

We now examine kinematical streaming motions in the data, 
i.e. rotation. Although there are data availabl e in the liter- 
ature for the stellar kinematics of NGC 4494 l|Bender et alj 
1 1994 hereafter B+94), we have obtained new long-slit data 
as part of a programme to obtain extended stellar kinematics 
in ETGs. These data will be presented in detail in a forth- 




150 200 
R (arcsec) 

Figure 5. Radial surface density profiles of the field stars (V- 
band; blue star symbols) and of the PNe (green squares) in 
NGC 4494. The PN number counts have been corrected for spa- 
tial incompleteness, and arbitrarily normalized to match the stel- 
lar data. The vertical error bars of the PN data in this and in the 
following figures represent the 1 a uncertainties (based in this case 
on counting statistics and completeness correction uncertainties), 
while the horizontal error bars show 68% of the radial range of 
the PNe in each bin. The purple solid curve is the Sersic model 
fit to the stellar photometry (see main text). 



coming paper (|Coccato et al.l [20081 ). and we summarize the 
main features of the NGC 4494 data here. 

Deep spectra were obtained with the ESO Very Large 
Telescope in service mode [Programme ID 076.B-0788(A)] 
using the FORS2 spectrograph, equipped with the 1400V 
grism, through a l'.'0-wide slit. The spectra were taken 
along the major (PA=0°) and minor (PA=90°) axes of the 
galaxy, and were re-binned in radius to reach S/N ~ 30 
at R ~ 100" ~ 2R a . Kinematical stellar templates were 
chosen from the In do-U.S. Coude Feed Spectral Library 
l|Valdes et alj |2004| ) and convolved with a Gaussian func- 
tion to match the instrumental FWHM. Various moments 
of the velocity distribution were extracted using the Penal- 
ized Pixel-Fitting method of lCappellari fc Emselleml (|2004h . 
The fitted spectral range was ~ 4560-5860A. 

The results for the stellar kinematics are shown in 
Fig. [6] The rotation profile along the major axis shows the 
kinematically decoupled core inside 5", as previously ob- 
served b y B+94, and coinci dent with the photometric stel- 
lar disc ( Carollo et al.|[l997l ). Since we are interested mostly 
in the outer kinematics, this very central feature is of no 
concern for us. 

Before examining the PN kinematics, we check their 
self-consistent systemic velocity by taking the median of the 
individual PN velocities in the azimuthally complete region 



8 Napolitano et al. 



E 40 5 



B 40 





E 100 - 



50 100 150 

R (arcsec) 



so ioo 150 aoo 

R (arcsec) 



50 300 




150 

(arcsec) 



Figure 6. Observed kinematics of stars (filled star symbols) and PNe (open circles) in NGC 4494, with uncertainties shown by vertical 
error bars. The data from opposite sides of the galaxy have been folded into a single radial axis. The horizontal error bars show 68% of 
the radial range of each PNe bin. The left panels show the mean velocities relative to the systemic velocity, and the right panels show 
the dispersions about the mean. The top panels show the major axis profiles, and the bottom panels show the minor axis. 



of 40"-135". This value is 1344 km s~\ exactly the same 
as the NED velocity for NGC 4494. 

We then measure the rotation of the PN system by cal- 
culating the mean velocity in bins within a wedge of 20 de- 
grees opening angle, aligned with the directions used for 
the long slit spectroscopy (PA=0° for the major axis and 
PA=90° for the minor axis). From Fig. [6] (left panels) it 
is apparent that the rotational properties of the PNe and 
stars are generally consistent. The discrepancy at small radii 
(~ 50") on the minor axis may be caused by the finite width 
of the PN measurement wedge (i.e., picking up off-axis ro- 
tatio n) and by statistical flu ctuations in the velocity sample 
(see lNapolitano et alJboOll ). 

The constant major-axis rotation velocity of the stars 
outside the galaxy core of 60 km s" 1 is found using the PNe 
to extend to at least 4 _R C . The minor-axis rotation is gen- 
erally small but does suggest some finite rotation at large 
radii, i.e. a kinematic misalignment. This feature could be 
caused by a transition in the orbit stru cture or by projection 
effect s in a smooth triaxial system (e.g. Ivan den Bosch et all 
2008), and may be related to the isophote twisting between 
~ 75" and ~ 200" (see Fig. |4j). 

In order to make use of all the PN data rather than just 
those along the principal axes, we construct 2-D maps of the 
velocity distribution, as in D+07. Assuming at least triaxial 
symmetry, we produce higher-S/N maps by mirroring the 
PN data points in phase space from (x, y, v) to {—x, —y, —v). 



We then smooth the velocity field using a median filter and a 
Gaussian kernel of width a ~ 60" (~ 4.6 kpc) on a regularly 
spaced grid of 75" cell widths (see also D+07). The resulting 
mean velocity field is shown in Fig. [7] (left panel) . Again we 
see evidence of rotation along the major axis (PA ~ 0°) in 
the galaxy's central parts, twisting to another axis (~ 20°) 
outside ~ 1 Re. 

We next try an alternative measure of the PN rotation 
that is insensitive to azimuthal incompleteness. It works by 
fitting the velocities to a cosinusoidal curve: 

v{4>) — v x cos(0 - <J) ), (2) 

where <f> is t he position angle of the PNe relative to the ma- 
jor axis (cf. iNapolitano et al.ll200lT ). We report the best-fit 
parameters for various radial bins in Table [4] Using this 
technique, there is no significant evidence that the rotation 
is misaligned with the major axis (PA=0°) at any radius. 
Although there is a twist at intermediate radii (consistent 
with the smoothed velocity field in Fig. , this is not sta- 
tistically significant and we will consider the whole sample 
as consistent with PA=0°(as in the first row of Table [4}. 

Finally, we wish to quantify the dynamical significance 
of rotation within 1 R c . This is best done with high- 
quality integral field st ellar kinematics data (e.g. SAURON; 
ICappellari et al.1l2007i . hereafter C+07), but we can make 
an estimate based on the available long-slit data. C+07 de- 
rived an empirical calibration (their eqn. 23) between the 
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Figure 7. Two-dimensional kinematics of PNe in NGC 4494. Individual PN data points are shown as X symbols, with values given by 
colors from the scale below. The data set has been doubled by a reflection symmetry, and then smoothed (see main text). Solid curves 
illustrate iso-velocity contours. The dashed ovals show the 2 R e isophote. Left: Mean velocity. Right: Velocity dispersion about the mean 
(the squarish contour is partially an artefact of the gridded smoothing technique) . 



Table 4. Best-fit parameters of sinusoidal rotation model to PN 
velocity data, in radial bins. <j>0 is the direction of maximum ve- 
locity. 



Radius ["] 


N. PNe 


vo [km s 1 ] 


00 [degrees] 


- 400 


255 


-23 ± 9 


-2± 25 


0-75 


98 


-28 ± 18 


24 ±39 


75 - 133 


79 


-12 ± 16 


42 ± 76 


133 -400 


78 


-37 ±14 


-8 + 19 



SAURON rotation dominance parameter within 1 R c 
the traditional measure of B+94: 



B+94 



and 



(3) 



where D max is the maximum observed rotation velocity along 
the major axis, and <jq is the average velocity dispersion 
within R e /2. Given v max = 76 km s _1 and Co = 163 km s _1 
from B+94, we find that (v/a) c ~ 0.26. Comparing this 
value to the isophote ellipticity from B+94 (e = 0.16), we 
find that NGC 4494 is located firmly in the space occupied 
by "fast rotators" (C+07 Fig. 11). Although the calibration 
above was done for the B+94 data only, very similar results 
would be obtained using our new photometric and long-slit 
data. 



3.3 Velocity Dispersion Profile 

We next examine the random motions in NGC 4494, repre- 
sented by the velocity dispersion profile. Fig. [6] shows that 



the PN and stellar dispersions are consistent along the ma- 
jor and minor axes. The dispersion declines steadily with ra- 
dius, and random motions dominate streaming motions out 
to ~ 100" (~ 2 R e ). Fig. (right) shows that the dispersion 
is roughly constant along the galaxy's isophotes. 

We reduce these data to a single velocity dispersion pro- 
file as a function of the intermediate radius R m . The rotation 
and true dispersion profile are folded into an root-mean- 
square velocity profile wrms = s/v 2 + a 1 , where v and a are 
the rotation and dispersion components respectively. The 
RMS velocity is a measure of the total kinetic energy, and 
we henceforth loosely refer to it as the velocity dispersion or 
VD. We combine the stellar data from the different axes by 
averaging, while folding the (small) systematic differences 
into the final uncertainties 

The resulting dispersion data are plotted in Fig. [8] The 
PN dispersions are similar to those reported in R+03, but 
are based on 3.5 times as many velocities, and extend to 
twice the radius. In our subsequent dynamical models, the 
stellar dispersion will be fitted out to 60" (~ 1-R C ), where 
it is more accurate, and the PN dispersion outside this ra- 



5 The uncertainties in the PN dispersion use a classical ana- 
lytic formula that assumes a Gaussian distribution, i.e. Aiirms ~ 
yJHiV^ /2N . Howev er, we expect this to pro duce accurate results 
in realistic systems llNapolitano et aTfe oOl), and we have carried 
out additional Monte Carlo simulations of a simplified galaxy with 
radial orbits, finding that the dispersion is very accurately recov- 
ered with our estimator, with a possible bias to be ~ 5% too high. 
Also, the expression for stellar vrms 1s valid for comparison with 
the azimuthally-scattered PNe only because we are averaging the 
major and minor axis long-slit data; if only major axis data were 
available, the expression would be fRMS = -\/f 2 /2 + cr 2 . 
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Figure 8. Composite projected velocity dispersion profile of 
NGC 4494, with data from stars (filled star symbols) and PNe 
(open circles). The solid curve shows a power-law fit to the PN 
data outside 60". 



dius. To characterize the decline of the dispersion with ra- 
dius, we fit a power-law to the PN data outside 60", using 
a maximum-likelihood method that fits the discrete data 
without binning, though it does assume an unde rlying Gaus- 
sian velocity distribution l|Bergond et al.l [20061 ). We find a 
power- law exponent of —0.2 ±0.1 (see Fig. [SJ, which is sig- 
nificantly shallower than -0.6 ± 0.2 found for NGC 3379 
(D+07), suggesting a dynamical difference between the two 
galaxies (discussed further in Q. 

We also test the sensitivity of the derived velocity dis- 
persion profile to the outliers selection (Sec. 12.2] ). Restoring 
all the objects to the PN sample which were rejected at 
n = 2.5 but not n — 5, we find that the only affected feature 
is a peak in the dispersion between 100" and 150", which 
does not change the broad trend. 

3.4 Higher-Order Velocity Moments 

The velocity dispersion is a crude measure of the particle 
speeds in a collisionless system, whose finer orbital proper- 
ties can be probed using higher-order moments of the line- 
of-sight velocity distribution (LOSVD). Stars on isotropic 
orbits in a gala xy halo with a flat rotation curve have a Gaus- 
sian LOSVD (|Gerhard| [T993). In a tangentially-anisotropic 
system where the particles follow nearly circular orbits, the 
observed velocities tend to pile up at the circular speed, 
producing an LOSVD which is flatter than a Gaussian, or 
even double-peaked. Radial orbits produce centrally peaked 
LOSVDs with long tails. In the case of integrated-light stel- 
lar velocities with fairly high S/N, these LOSVD shapes 
are described via the Gauss-Hermite momen ts such as /14 
l|van der Marel fc Franxlll993l ; lGerhard|[l993l '). In the case 
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Figure 9. Projected kurtosis profile of NGC 4494, with symbols 
as in Fig. [8] 



of discrete velocities, it is more straightforward to compute 
the classical dimen sionless kurtosis, k = v 4 /(v 2 ) 2 — 3 (see 
Ijoanes fc Gilllll~998l for exact expression and uncertainties . 
Broadly speaking, k ~ for isotropic orbits, « < for tan- 
gential orbits, and /-c > for radial orbits. 

We quantify the shapes of the stellar and PN LOSVDs 
in NGC 4494 in Fig. [ 51 using the approximate re lation 
()van der Marel fc Franxll 19931 ). The 



SVEIh for the stars 



PN kurtosis is fairly noisy, but is consistent with the stel- 
lar properties in the region of overlap. However, there is a 
systematic offset between our stellar /14 values and those re- 
ported by B+94 inside 20", where we believe our improved 
template library procedures yield more reliable results. In 
the halo, the PN kurtosis is consistent with zero. 



4 DYNAMICAL MODELS 

We now combine the photometric and kinematical data for 
the stars and PNe in NGC 4494 into integrated dynamical 
models in order to derive the mass profile of the galaxy, and 
test whether or not it hosts a DM halo as expected from 
ACDM. In 94.11 we begin with the simple pseudo-inversion 
mass modelJ used in R+03 and D+07, and introduce a new 
series of physically-motivated kurtosis-based models in 94.21 
We address systematic issues in 94.31 

In both of our main modelling exercises ( 94. H and 94. 2|) . 
we will apply the spherical Jeans equations to NGC 4494, 



6 Again, the uncertainty calculation assumes a Gaussian distri- 
bution, while for a more general distribution the errors are for- 
mally infinite. However, the Monte Carlo simulations mentioned 
in EI3.3I demonstrate accurate recovery of the kurtosis using our 
estimator, with a systematic deviation of no more than ~ 0.1. 

7 This is a phenomenological mass invertion method where we re- 
construct the mass distribution from the velocity dispersion data. 
As shown in i]4.1l since we make use of a parametrized model for 
the o r and the anisotropy parameter in the radial Jeans equation, 
this is not a full inversion method. 
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Figure 10. Model fit to the NGC 4494 velocity dispersion data 
using the pseudo-inversion mass model and assuming isotropy. 
The best fit and 1 a range of uncertainty are shown by the solid 
curve and shaded region, respectively. 



since the observed stellar isophotes are very round (see £|3 . 1 1 
for photometry, and D+07 §8.1 for discussion of the spherical 
approximation). We will examine a non-spherical model in 
i]4.3l and Appendix [C] Given the consistent agreement seen 
in §S] between the stellar and PN properties, we will assume 
that all of these data are drawn from the same underlying 
dynamical tracer population. We will also in general omit 
the stellar kinematics data inside 10" from our model fits, 
since there appears to be a strong dynamical change in this 
nuclear disc region which our smooth Jeans models are not 
designed to reproduce. 

4.1 Pseudo-inversion mass model 

Our first series of models uses simple assumptions to convert 
the observed kinematics into a mass profile M(r). This phe- 
nomenological approach, as introduced in R+03 and D+07, 
is computationally light, need not involve Abel inversion 
integrals, and has no prejudice about the form that M(r) 
should take, nor about the stellar M/L value (which will 
be discussed later in this section). On the other hand, it is 
not possible to directly test any theoretical predictions for 
the DM distribution, and the resulting mass profile may not 
even be physical. 

The mass modelling procedure consists of five steps: 

(i) Adopt a simple smooth parametric function for the 
intrinsic radial velocity dispersion profile: 



oy(r) = (T 



1 + 



(4) 



(ii) Assume a fixed anisotropy profile: 
(3{r) = 1- al/al, 



(•») 



where ag and oy are the tangential and radial components 
of the velocity dispersion ellipsoid, expressed in spherical 
coordinate^]. 

(iii) Project the radial component of the spatial velocity 
dispersion ay (see Eq. IB4[I for comparison with the line-of- 
sight velocity dispersion data <ji os (R). 

(iv) Iteratively adjust the free parameters in equation |4j 
to best fit the model to the observed dispersion profile. 

(v) Use the best-fit model (Eq. [4j in the Jeans equation 
4-55 of lBinnev fc Tremainel l| 19871) to calculate M(r): 



M(r) 



G 



d In j, d In ay 
dlnr 



dlnr 



+ 2/3 



(6) 



where ao,ro,rj are a minimalistic set of free parameters. 



where j*(r) is the spatial density of the PNe, and corre- 
sponds to an Abel deprojection of a smoothed density law 
fitted to the stellar data. Additional quantities may then be 
computed, such as the cumulative M/L. 

Starting with the minimalist assumption of an isotropic 
galaxy (/3 = 0), we find that the simple model @ is able to 
fit the dispersion data very well (Fig. llOp . The resulting M/L 
profile increases clearly with the radius (Fig. Hip , providing 
a strong indication for the presence of an extended DM halo. 
Using 5J?o (241") as a benchmark value, we find that the 
B-band M/L within this radius is Ts,s = 6.3lg® T©^, 
where quoted errors account for the 1-a confidence region in 
the parameter space (ao,ro,rj) of the dynamical model. As 
we will see in t|4.2l the derived value for T5 is sensitive to 
assumptions about /3(r) but not about the detailed form for 
M(r), and so is useful as a robust quantity for comparison 
to theory. 

The shaded regions in Figs. llOl and llll show the statisti- 
cal modelling uncertainty, but there remains the systematic 
uncertainty implied by the anisotropy in step (ii) of the pro- 
cedure. Adopting constant values of (3 = ± 0.5 for a plausible 
(though not exhaustive) range of the stellar anisotropy, we 
can find fits to the data just as good as in Fig. 1101 However, 
the famous mass-anisotropy degeneracy implies very differ- 
ent M(r) profiles, as shown in Fig. 1111 Assuming (3 — +0.5 
implies more DM than with (3 = 0, while (3 = —0.5 implies 
a smaller but non-zero amount of DM (i.e. a constant M/L 
is excluded at the 1 a level). Note though that the kurtosis 
values found in i)3.4l suggest (3 <; and thus that the higher- 
mass solutions are to be preferred (as will be explored in the 
next section). 



8 Since we are folding the observed rotational and pressure sup- 
port into a single parameter drmSi we use a similar approach 
in the Jeans equations, so that technically /3(r) si— (cr^ + 

irp/(2<r;?). We thus assume that Vg = ?;? everywhere, in contrast 
to an alternative assumption t hat = cr^, wi th the rotation sup- 
port modelled separately (e.g. Irlui et a,l.|[l995h . Strictly speaking, 
a spherical model is not warranted for a rotating system any- 
way, but in the limit of weak rotation, it is a matter of adopting a 
convenient dynamical approximation. The rotational contribution 
does rise in the outer parts of NGC 4494 to as high as v/a ~ 1.4, 
potentially causing systematic issues in our mass analysis which 
would require a more generalized study to quantify. 
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Figure 11. Cumulative B-band mass-to-light ratio (M/L) of 
NGC 4494. Most of the curves show results from the mass- 
inversion method: the black solid curve with grey shaded region 
shows the best-fit (3 = model with its statistical uncertainties; 
the orange dot-dashed and blue dotted curves show the best fit 
models with (3 = —0.5 and /3 = +0.5, respectively; the green short 
dashed curve shows the adopted /3(r) model of Eq. [7] The purple 
dot-long-dashed line shows the best cosmologically-based model 
["NFW + /3(f)"] and the red double-dot-dashed curve is the log- 
arithmic potential model with /3(r) (see £14.2.51 and 34.2,61 1. The 
horizontal blue long-dashed line with surrounding shaded region 
shows the stellar M/L and its uncertainty. 



As a final refinement, we wish to test an anisotropy pro- 
file based on theoretical expectations. In general, collision- 
less systems such as galaxy haloes and clusters that were 
formed through mergers should have some degree of radial 
aniso t ropy (e.g..lDiaferidll999l;IColm et al.ll200ol ; lRa 



et al.l 

120041 ; IWoitak et al.ll2005l : lAbadi et al] I2006T I— a prediction 
that has found empirica l support from the Galactic halo and 
from elliptical galaxies (|Chiba fc B eers 200fj; iGerhard et al.l 
l200ll ; DL+08a,b). In particular, ML05 have introduced a 
radially- varying anisotropy profile, 



0(r) = A) 



r + r a 



(T) 



finding /3o — 0.5 and r a ~ 1.4 .Re when applied to the merger 
simulations of D+05. Adopting this profile with r a = 65", 
we find a M/L profile that is slightly shallower than the 
pure f3 = +0.5 model (Fig, in). 

The overall plausible range for the benchmark-radius 
M/L of NGC 4494 is T 5 ,s = 5.0-9.1 T Q , B (including both 
statistical and assumed systematic /3 uncertainties). This re- 
sult is very similar to the pseudo-inversion mass approach 
from R+03, which was based on a different photometric 
model and a much more limited data set: Ts.s = 5.6- 
8.0 Tq,b (distance assumption revised, and only the f3 un- 
certainty included). 



To further contextualise our result, we calculate the 
M/L gradient parameter introduced by N+05: 



R C AT 
T in AR' 



(8) 



where Ti n is the central dynamical M/L. This parameter 
is independent of the bandpass and the distance, and is 
fairly insensitive to the modelling details within the spheri- 
cal Jeans formalism (with the exception of anisotropy) . For 
simple interpretations we can assume that Tin is dominated 
by the stars as in the case of NGC 4494, but independently of 
this assumption, the VfY parameter has been found to show 
strong observational correlations with other galaxy proper- 
ties (N+05). For NGC 4494 we find V^T = 0.02-0.23. This 
result is nearly identical to our previous finding, placing the 
galaxy among the po pulation of low-V^T (i.e. low-DM) sys- 
tems (see N+05 and lFerreras et al.ll2005l ). 

For NGC 3379, the mass-inversion results of D+07 im- 
ply V^T = -0.03-0.14. The two galaxies thus appear to 
have consistent amounts of DM given the total uncertain- 
ties, but assuming the same anisotropy profile for both, the 
two galaxies are just barely consistent at the 1 a level. The 
reason for this difference is that the dispersion profile of 
NGC 4494 is shallower than that of NGC 3379 (see 333)) . 
while its stellar luminosity falls off more steeply (which is a 
factor that should steepen the dispersion profile). 

We can also independently estimate the amount of 
mass contained in the stars (i.e. the stellar M/L, T*) 
via independent modelling of the stellar population based 
on line strengths. Unfortunately, there is not a radially- 
extended population map of NGC 4494 available as for 
NGC 3379, and we are forced to resort to constraints on 
its very central regions. Based o n luminosity-weigh ted line 
strengths r eported in s ide ~ 1.2" (|Trager et al.lfl998l ). we use 
the on- line IWorthevI i| 19941 ) modelfl to estimate an age of 
~ 10.6 Gyr and a metallicity of [Fe/H] ~ +0.21, implying 
T,, s ^ 4.3 ± 0.7 T , s (Kroup a IMF assumed, ari d ~ 20% 
typical age errors adopted from lTrager et al.llBo O'Fl. while 
higher values would be obtained using a Salpeter IMF. A 
caveat here is that central line indices in general are very 
likely to be contaminated by young stellar subpopulations 
confined to the central regions (e.g. in the nuclear disc of 
NGC 4494) and thus to produce lower T, values than in the 
halo regions that concern us. 

We plot the T* value in Fig. 1111 and see that it is simi- 
lar to the dynamically-inferred M/L in the galaxy's centre. 
Although there are large uncertainties in this comparison, 
the coincidence of the best-guess M/L values suggests that 
the stars contribute the large majority of the mass in the 
central regions (similar conclusi ons were found in th e case 
of NGC 3379 by D+07; see also IGerhard et al]|200l] ). It is 
only outside ~ 2 R e where the effects of DM become no- 
ticeable. In the next section we will attempt to constrain T* 
by dynamical means, and to reconstruct the DM profile in 
more detail. 



9 http:/ /astro. wsu.edu/worthey/dial/dial_a_model. html 

10 Somewhat l ower T « values are implied by t he results of 
IGerhard et al] ll200lh . iDenicolo et al] (20051) and IZhang et all 
J2008h 
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4.2 Multi-component models 

We now move beyond the pseudo-inversion mass method 
used in R+03, D+07 and 8 14.11 and construct physically- 
based dynamical models of NGC 4494, allowing for more 
direct comparisons of observation and theory. Along with 
these models, we include a kurtosis-based application for 
alleviating the mass-anisotropy degeneracy. This approach 
has similarities to the orbit models used for NGC 3379 in 
R+03, as we will discuss in £14.2.31 It is outside the scope of 
this paper to construct fully rigorous and general dynamical 
models. Rather, our aim is to generate plausible "best-guess" 
models that reproduce the broad features of the data. 

Our two-component mass model is simple but realis- 
tic, consisting of a luminous field star distribution and (op- 
tionally) a DM halo. The total gravitational potential may 
thus be expressed as $ = $* + 4>d- The stellar gravitational 
potential <J>*(r) is derived from the stellar luminosity j»(r) 
(according to the Sersic model of £|3.1[) . combined with some 
assumed constant T,. 

We first describe our suite of mass models, where the 
DM distribution follows either the NFW profile (" 34.2. l|l or 
a pseudo- isothermal form f £|4.2.2[) . Our dynamical methods 
are outlined in £14.2.31 and our results given in £)4.2.4tf4.2.6l 



4.2.1 NFW model 

Our first suite of mass models is based on simulations of 
collisionless DM halo formation in a ACDM cosmology. In 
this case the DM density takes the approximate form of an 
NFW profile: 



Pd(r) = 



{r/rs){l + r/r a y , 



(9) 



where p s and r s are the characteristic density and scale ra- 
dius of the halo. The cumulative dark halo mass is 



M d (r) = 4np a rlA(r/ra) 
where 

A(x) = ln(l + as) - 
The potential is: 



1 + x 



*d(r) = 



4yrGp s 



In 



(—) 

\r + r s J 



(10) 



(11) 



(12) 



where G is the gravitational constant. 

The three free parameters describing the NFW mass 
model are thus Y*, p s and r s . The halo can alternatively be 
parametrized by the virial mass and concentration, M v i r = 
47rA v i r p cr i t r5i r /3 and c v ir = r v i T /r s , where the critical den- 
sity is p cr it = 1.37 x 10 _ Mq pc -3 and the virial overdensity 
value is A v i r = 101. The expected values for these model pa- 
rameters are not arbitrary in ACDM — a theme to which we 
will return in SjS] 



4.2.2 LOG model 

Our second mass model co nsists of a logarithmic potential 
|Binnev fc Tremainel Il987l §2.2.2) which was motivated by 
observations of spiral galaxy rotation curves. The potential 
is: 



(13) 



where vo and ro are the asymptotic circular velocity and 
core radius of the halo. The corresponding DM density and 
cumulative mass profiles are respectively: 



Pd(r) 



and 



M d (r) 



^o 2 (3rg+r 2 ) 
4%G{rl + r 2 ) 2 



v 2 r 3 



G r 2 + r 2 



(14) 



(15) 



The three free parameters of this "LOG" model are thus 
T„, vo and ro. We define a virial mass according to the 
same definition as in £14.2,11 

Unlike the cuspy r _1 density core of the NFW halo, 
the DM density of the LOG halo is constant in the centre. 
Outside of the core, the density decreases as r~ 2 to produce 
a constant circular velocity profile v c = GM(r)/r , similar 
to the NFW model near r = r s . Such a model allows us to 
maximize the stellar contribution to the mass in the cen- 
tral regions, and thus test a "minimal DM halo" scenario. 
Similar models ha ve been used to successfully fit galax- 



1991 




2007 





after DL+08a; DL+08b). 



4-2.3 Dynamical methods 

The Jeans modelling approach employed here is more tradi- 
tional than that of £14.11 one starts with a trial mass model, 
solves equilibrium equations for the internal velocity mo- 
ments, and projects those moments onto the sky for com- 
parison with the data. This approach typically means solv- 
ing for the second-order velocity moments (i.e., the veloc- 
ity dispersions) and comparing to projected velocity disper- 
sion data. However, such data cannot constrain the mass 
model without additional strong assumptions; the entire 
LOSVD in a spherical system is necessary to have any hope 
of uniquely determining its phase-space distribution (e.g. 
iDeionghe fc Merritj|l992i ). In practice, the systematic un- 
certainties can be strongly reduced by using fourth-order 
moments of th e LOSVD to diagnose orbit anisotropies (e.g., 
lGerhard|[l993l ). 

One such method is to construct fourth-order Jeans 
equations in addition to the usual second-order, and to con- 
strain them with kurtosis data (e.g. Magorrian fc Ballantvnei 
l200ll : lLokasI 120021 ; lLokas fc Mamonl 120031 "). Although the 
higher-order Jeans equations are not closed in general, one 
can adopt a simple choice for the distribution function which 
makes the problem tractable (see Appendix|B}. This simpli- 
fication is arbitrary (e.g. f3 is assumed to be constant with 
radius) and does restrict the generality of our results, but the 
model is still more general than an assumption of isotropy. 
Given these restrictions, we can now constrain /3 using the 
data. 

The basic steps of our analysis are as follows, with de- 
tails for steps (ii) and (iii) provided in Appendix [B] 

(i) Set up a multi-dimensional grid of model parameter 
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space to explore, including j3 and the mass profile parame- 
ters (T*,p s ,r s ) or (T t ,v ,r ). 

(ii) For each model grid-point, solve the 2 nd - and 4 th - 
order Jeans equations. 

(iii) Project the internal velocity moments to <ti os and 

(iv) Compute the \ 2 statistics, defined as 

J Vdat a r oba mod "| 2 

x = E K^M ■ ^ 

where p° hB are the observed data points (<ti os and ki os ), 
pf 10 * 1 the model values, and Sp° hB the uncertainties on the 
observed values, all at the radial position Ri. We fit the 
PN data outside 40" and the stellar dispersion between 20" 
and 60". The stellar kurtosis is used between 30" and 60". 

(v) Find the best fit parameters minimising the x 2 ■ ^ n 
practice, we find that the VD is affected by both the mass 
and anisotropy profiles, while the kurtosis is almost entirely 
driven by the anisotropy. 

For comparison, R+03 modelled NGC 3379 with a sim- 
ilar suite of spherical mass profiles in non-parametric orbit 
models (A.K.A. "Schwarschild's method") that included di- 
rect fits to the stellar hi moments as well as to the full PN 
LOSVDs. Such methods (see also lChaname et alj|2008l and 
DL+08b) may be more powerful than our Jeans approach, 
but the latter is computationally faster and somewhat more 
intuitive. 

4.2.4 Results: no-DM case 

We start with the simplest model, assuming a pure-stellar 
potential (p 3 = or vq — 0) and isotropy for the veloc- 
ity ellipsoid (/3 = 0). The best-fit parameters of the model 
are given in Table \S\ together with their \ 2 values. Given 
the freedom to adjust T„, this model can fit the VD in the 
central regions ( <, 2R C ), but falls off too quickly in the 
outer regions (Fig. If 21 upper left). Although the kurtosis is 
included in the fit (Fig. 1121 upper right) , omitting these con- 
straints would not provide an appreciably better fit to the 
VD, since the M(r) and /3 assumptions already specify fully 
the shape of the VD. Note that in R+03 we were able to 
fit a pure-stellar model to the VD data, with the difference 
caused by a slightly steeper new stellar luminosity profile, by 
a ~ f 0% lower VD normalisation from the new stellar kine- 
matics data, and by a shallower outer VD slope as derived 
from our larger, more extended new PN data set. 

In order to increase the VD at large radius, we exper- 
iment with introducing a degree of tangential anisotropy 
(ft < 0). As seen in Fig. If 21 this strategy does improve the 
VD fit (best-fit f3 = — f.3), but requires a negative kurto- 
sis in the halo that is ruled out by the PN data. Given the 
observational uncertainties in the stellar luminosity profile 
j*(r) ( §3.1 [) . we test out a shallower outer profile (Sersic in- 
dex m — 5.2j"l. We find that we can indeed improve the 

11 Our surface photometry is originally g-band transformed to 



fit to the VD, but the kurtosis is still problematic for this 
model. The use of higher-order velocity moments therefore 
leads us to conclude that there is a non-zero amount of DM 
in NGC 4494 (as also found by R+03, D+07 and DL+08b 
for NGC 3379). 

Given the presence of DM in this system, we next 
determine what halo parameters are best consistent with 
the data for the two assumed DM profiles. 



4.2.5 Results: NFW model 

Once again we start by assuming isotropy, and find a best fit 
as shown in the central panels of Fig. If 21 with parameters 
again reported in Table [5] Note that inside ~ 50" ~ R c , 
both the stars-only and the stars+DM models can recover 
the VD profile with simple assumptions, suggesting that in 
these regions it is the stellar mass that determines the main 
kinematical features. 

The isotropic NFW solution is a reasonably good match 
to the data, but it can be seen that the predicted outer VD 
is somewhat too high, and the predicted kurtosis somewhat 
lowB 

Therefore we next allow for the anisotropy (3 to be a 
free parameter, though constant with the radius. The best-fit 
solution has a \ 2 value decreased by 9, although there is only 
one less degree of freedom. This implies that the anisotropic 
model is strongly preferred over the isotropic one, as also 
visible in Fig. If 21 where both the VD and the kurtosis are 
at last reproduced well at all fitted radii. The solution has 
a moderate degree of radial anisotropy (/3 ~ 0.5), and is 
formally preferred at the 3 a level over the best anisotropic 
stars-only solution; the inferred mass profiles are dramat- 
ically different, with the NFW model v c remaining much 
flatter with radius than the stellar model (see Fig. If 3[>. The 
radial anisotropy decreases the VD at ~ R c and requires a 
slightly higher T* value than in the isotropic case. 

The halo parameters in this anisotropic analysis are il- 
lustrated by Fig. If 41 which shows the joint region of per- 
mitted values for r s and p s , marginalized over the other free 
parameters, T* and (3o- It can be seen that there is a lim- 
ited degeneracy between the parameters, with alternatively 
low-p s and high-r s models (corresponding to high M v i r and 
low c v ir), or high-p s and low-r^, (low M v - lr and high c v ir). 

The best-fit anisotropy value in our modelling of /3 = 
0.46 is remarkably similar to theoretical expectations for 
galaxy haloes (see fl4.ip . Therefore in order to make an im- 
proved estimate of a plausible ACDM halo that would re- 
produce our data, we adopt Eq. [7] for the anisotropy pro- 
file /3(r). This model is isotropic in the centre and radially 
anistropic in the outer parts, where the transition radius 
r a is free to vary (within the range 30" to 110"). Because 
our Jeans modelling procedures are not entirely designed to 
cope with a radially- varying (3 (see Appendix [B}, we adopt 
a slightly different procedure than the preceding. We fit 
the VD data as before, but fit the kurtosis data only at 



the B- and V-bands, so the true stellar mass density profile (~ K- 
band) is likely to be steeper. 

12 Although a galaxy halo with (3 = and constant v c should 
result in Ki os = 0, our mass model deviates enough from this 
picture to produce non-zero kurtosis. 
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Figure 12. Multi-component Jeans model fits to the NGC 4494 kinematics data. The stellar data are shown by star symbols, while the 
PN data are open circles. The left panels show the projected RMS velocity profiles, while the right panels show the projected kurtosis. 
The top panels show fits from self-consistent stars-only models, the central panels show stars+NFW halo models and the bottom panels 
show stars+LOG halo models. The curves correspond to models as in the panel legends. See text for details. 
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Table 5. Summary of best-fit multi-component model parameters. 



Model 


ft 1 


( T J,-b) 


log M, 3 
(M @ ) 




log i\/ vir 5 
(Mq) 


/v, r 6 


/dM,5 7 


T( Rc) 8 

( r e,B) 


r B5 9 

( T 0.b) 


Y(-R vir ) 10 
( T 0,b) 




xVd.o.f. 12 


No-DM model 


star iso 
star tang 



-1.3 


4.6 
4.4 


11.04±0.06 
11.02±0.06 




11.04±0.06 
11.02±0.06 










4.6 
4.4 


4.6 
4.4 


4.6 
4.4 




(i 


25/37 
27/36 


NPW model 


NFW iso 
NFW+/3 
NFW+/3(r) 




0.46±0.15 
0.43±0.05 


4.1 
4.3 
4.1 


10 99+ ' 05 
iuMM -0.07 

11 01+ ' 05 
i± u -0.06 

10 99+ 05 
-0.0T 


6 9+ 05 
D -0.4 

6 5+ 1 ' 6 

s , +0.9 

8 - 3 -0.7 


12.03±0.12 

12 17+ 016 
-0.20 

12.05±0.10 


10±3 
13±5 
11±3 


0.41±0.08 
0.42±0.07 
0.44±0.08 


4 5+°' 2 
-0.1 

4 7+0-3 
' -0.2 
4 6+0-2 
-0.2 


«±; 


45+1° 


0.14 
0.15 
0.14 


21/35 
10/34 
5/22 


LOG model 


Model 


ft 1 


(T ,b) 


log M, 3 
(%) 


V 13 
1 
(kms -1 ) 


log M vir 5 

(Mq) ( 


arcscc) 


/dM ; 5 7 


T(ii c ) S 

CT ,s) 


r B5 9 
(T ,b) 


T(fl vir ) 10 
(T ,B) 


V^T 11 


x 2 /d.o.f. 12 


LOG iso 
LOG+ft 
LOG+/3(r) 




0.44±0.10 
0.43±0.07 


4.3 
4.3 
4.3 


11 01+ 007 
11 01+ ' 05 
11 01+ ' 05 


147 
172 
146 


12.19±0.18 
12.39±0.15 
12.19±0.18 


244 
205 
190 


0.34±0.14 
0.45±0.13 
0.39±0.16 


4 4+0-1 
-0.1 
4 4+0-2 

4 4+0-2 
-0.1 




^ifo 

100+|° 


0.11 
0.18 
0.15 


12/35 
10/34 
4/22 



NOTES: 1) Anisotropy at the benchmark radius of 5 R c ; 2) dynamical stellar mass-to-light ratio M/L, in S-band Solar units: typical uncertainty 
is ±0.2Tq ^; 3) log of stellar mass in solar units; 4) concentration parameter (see fl4.2. 1\ ; 5) log of virial mass; 6) ratio of total dark and luminous 
matter within the virial radius, f v ^ T — M^/M* at T" vir ; 7) dark matter fraction, /dm — -^d/C^d + a * 5 J? c ; 8) dynamical M/L at R e ; 

9) dynamical M/L at 5 R e ; 10) dynamical M/L at the virial radius; 11) M/L logarithmic gradient; 12) x 2 statistic (see text for details of data 
included); 13) asymptotic circular velocity: typical uncertainty is ±30 km s — 14) halo core radius: typical uncertainty is ±50 // . 
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Figure 13. Radial mass distribution of NGC 4494. The left panel shows the cumulative mass, and the right panel shows the circular 
velocity profile. The heavy solid curves are the best-fit stars+DM model with varying anisotropy (black is NFW, gray LOG) while the 
light black solid line is the isotropic NFW case; the dashed curve is the stellar distribution for the LOG model (Y„=4.3Yq g). The green 
dotted curves show the fi = 0.5 stars+NFW model and the red dot-long-dashed curves the (3 = 0.5 stars+LOG model while the blue 
dot-dashed curves show the (3 = 0.5 pseudo-inversion mass result. A legend for the model symbols is in the left panel. 



large radii (R > 100"), where the anisotropy is fairly con- 
stant, and considered fixed to its value at 200" (~ 47? e ); see 
eqn. <[B§]|, 

This final best-fit theoretically-motivated model is 
again a good match to the data (see Fig. [12] and \ 2 m Ta- 



ble[5]). The anisotropy transition radius is r a = 45"±20", i.e. 
radial anisotropy may set in at smaller radii in NGC 4494 
than in the simulations of D+05. The stellar M/L of this 
solution of T»,b = 4.1 ± 0.2 Tq,b agrees remarkably well 
with the independent estimate from stellar populations of 
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Figure 14. Halo parameter solutions for NGC 4494 NFW-based 
model, plotted as halo scale radius r s vs central density p s . The 
contours show the two-parameter \ 2 confidence levels (1 and 2cr). 
Dashed contours show the constant-/3 model, and solid contours 
the /3(r) model. The best-fit solutions are marked with crosses. 
The solid line with shaded region shows the theoretical expecta- 
tions for the halo parameters (see JJeOJi and the 1— a scatter. 

T„, B ~ 4.3 Tq, s ( 34~Tj) . The central NFW halo parameters 
of p a = 0.0019 ± 0.0006 M Q pc" 3 and r s = 410" ± 30" = 
32 ± 2 kpc correspond to a virial radius, mass and concen- 
tration of r vir = 261 ± 30 kpc, M vil = (1.12t°H) X 1O 12 M 
and Cvir ~ 8 ± 1. 

As shown in Fig 1141 this f3(r) solution is better con- 
strained than the constant-/? one, and coincides with its 
secondary \ 2 minimum — showing that the results were not 
sensitive to the potentially worrisome stellar kurtosis points 
between 30" and 60". The data in both models are incon- 
sistent at the ~ 1 a level with the predicted relation from 
ACDM (shown in the same Figure, see Eq. [19] in 35. 2[) . As 
a visual demonstration of the significance of this result, in 
Fig.[T2]we show also a "typical" expected ACDM halo, with 
c vir = 12.3 and M vir = 2 x 10 12 M Q , and a f3(r) profile with 
r a = 67". This model is a relatively poor fit to the disper- 
sion profile, and required a smaller T»(=3.4Tq i b) in order 
to match the central data points. On the other hand, the 
kurtosis is well matched, as the anisotropy profile is similar 
to our best-fit NFW+beta(r) model. 

4-2.6 Results: LOG model 

We next carry out an equivalent model sequence for the LOG 
halo model, with results shown in Fig. [12] (bottom) and Ta- 
ble[5] For the cases of (5 — or (5 =const, the LOG model can 
fit the data significantly better than the NFW model, owing 
to the former's slower transition with radius from stellar to 
DM dominance. In the case of variable anisotropy, the LOG 
and NFW models have comparably good fits and similar 
/3(r) profiles, but in the context of the LOG model, we can- 



not distinguish between the different /3 models. However, we 
will adopt the /3(r) solution as the best LOG model because 
of its theoretical motivation. 

The M/L at 5 R c is very similar for the f3(r) NFW and 
LOG models, demonstrating that this quantity is well con- 
strained by the data, independently of the details of the mass 
models. The mass profiles are also similar given the mod- 
elling uncertainties (Fig. ll3[) , with the NFW model assuming 
a relatively low concentration which mimics the constant- 
density core of the LOG model, and with the NFW scale 
radius value of r s ~ 9 R a producing a roughly flat v c profile 
in the halo regions as in a LOG model. The remaining differ- 
ences in the model profile shapes result in slightly different 
values for T», which is lower for the NFW model because of 
the stronger central DM contribution from its cuspy core. 

Although we cannot at this stage distinguish between 
the NFW and LOG models, there are possibilities for do- 
ing so in the future. New data extending to larger radii 
could verify our current results and provide stronger dis- 
criminants between the mass models surveyed here and al- 
ternative ones — but would not be able to discriminate be- 
tween the NFW and LOG models, which begin to diverge 
only near the virial radius (note the differing extrapolated 
M v i r values in Table [SJ. More crucial would be generalized 
modelling which fits the fine details of the stellar kinematics 
in order to ease out the differences in the total mass pro- 
file shapes in the central regions. In this context, additional 
data over a wide baseline in radius would also help to build 
up the S/N of the constraints. 

4.3 Systematics and mass comparisons 

We next investigate the sensitivity of our mass results to 
some systematic issues and to different modelling methods 
and mass tracers. We can first examine the results from the 
different modelling methods and assumptions we have used 
thus far. As seen in Figs. [TT] and 1131 the total mass profile 
at the outer extent of the data (~ 300") is fairly well de- 
termined (within ~ ±30%) regardless of the anisotropy and 
mass model assumptions, and of the Jeans method used (cf. 
also the consistency of the Tbs benchmark parameter in 
Table [S]). Similarly, the inferred value for T» is the same at 
the ~ ±5% level, and is reassuringly consistent with the in- 
dependent estimate from stellar populations analysis. Using 
the more detailed models of H4. 2 1 with their anisotropy con- 
straints, we update our estimate of the M/L gradient ( H4.1|l 
to V<T ~ 0.1-0.2. 

One caveat about all of these Jeans models is that a 
somewhat arbitrary radial region is chosen for fitting the 
central stellar kinematics (thereby normalising the entire so- 
lution), since the models are not sufficiently flexible to re- 
produce the high-5/A bumps and wiggles representing fine 
structure in the s tellar data. The only such model to-date of 
NGC 4494 was bv lKronawitter et af] (|2000h . based on spher- 
ical non-parametric distribution function fitting. Using the 
B+94 stellar kinematics data-set, they found a circular ve- 
locity profile at ~ 0.5 R a which is slightly different from ours 
(perhaps owing to the data differences: see £|3.4[1 . Similar re- 
sults without the benefi t of constraints were found by 
iMagorrian fc Ballantvne] (|200ll ). 

It is of critical importance for more fully understand- 
ing this system to employ a more general dynamical model 
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that can fit the stellar data in detail, while not requiring the 
Lokas-type simplification (see Appendix[B} we have imposed 
on the distribution function. There are already some indica- 
tions from our Jeans modelling that other regions of model 
space may be permitted (with lower T* and higher central 
DM content). While more general modelling is outside the 
scope of this paper, we can estimate the possible system- 
atic effects from the Jeans smoothing simplifications. In the 
generalized NMAGIC modelling of NGC 3379 (DL+08b), 
the halo anisotropy could not be constrained to better than 
A/3 ~ ±0.3. Taking this as a guide, we re-run our NFW- 
based Jeans models with /3o fixed at either +0.2 or +0.8, 
finding that the corresponding systematic uncertainty in c v ir 
is then ±1.8, and in logAf v ir is ±0.26. 

Our models make powerful use of the assumption that 
the PNe and field stars are drawn from the same underly- 
ing dynamical population. This premise was questioned by 
D±05, who posited that the PNe could reflect a recent burst 
of star formation with distinct dynamics from the observed 
field sta rs — a hypothesis with some support from obser - 
vations (jSambhus et all 120061 ; iForestell fc Gebhardd I200II . 
We discussed this issue in more detail in D±07 and 
ICoccato et~al] (|2008T ). and here note that NGC 4494, al- 
though a^easc^rahlecandidate for a young stellar subpopu- 
lation dDenicolo et alj [20051 ). in fact shows good agreement 
between its PN and stellar properties (see e.g. ij3.ip . 

Our ACDM mass models assume a fairly basic colli- 
sionless NFW form for the DM distribution. This profile 
is bound to be altered by interaction with the collisional 
baryons; although the amount and even the sign of the effect 
are not clear, the standard view is that the change can be 
appro ximated by adiabatic contraction (e.g. iGnedin et al.l 
120041 '). The effect of this contraction would be to draw more 
DM into the central regions, flattening the rotation curve, 
and thereby requiring a lower stellar M/L (possibly los- 
ing the consistency of our solution with the independent 
value) and a lower concentration. A more recent revision 
of the collisionl ess DM distribution (the Einasto model: 
iGao et all I2008T ) would tend to push things in the oppo- 
site direction, but more weakly. Of course, these conven- 
tional mass-model pictures could be completely wrong, as 
we have seen that a cored-halo mode l fits our data equally 
well (see also- e.g. Ide Blok et alj|200ll : iBorriello et ail 12003 ; 
iGentile et all 12007^ 7 and alternative gr avity models should 
also be explored (e.g. iTiret et al.ll2007n . 

A final concern, and a perennial problem in mod- 
elling any galactic system, is the symmetry assumption. 
The spherical symmetry adopted here has in particular been 
called repeatedly into question, and certainly more robust 
results should eventually be derived using less restrictive 
symmetries and more detailed modelling. NGC 3379 is a 
similar case to NGC 4494, where its round apparent shape 
could be produced by a very flattened system seen near 
face-on. DL+08b used flattened particle-based "NMAGIC" 
models of NGC 3379 to explore the simultaneous effects of 
shape and anisotropy on the mass inferences from PNe. In- 
triguingly, they found that relaxing the spherical assumption 
hardly affects the inferred halo mass: the intrinsic flattening 
and inclination are strongly limited by the observed rota- 
tion field, and the main source of systematic uncertainty 
remains the halo anisotropy (which the available data do 
not constrain). Thus the effects of asphericity on mass mod- 



elling may not be a s important as previously thought (cf. 
iBekki fc Pen3l2006l ). although the effects of triaxiality have 
yet to be fully explored. 

It is beyond the scope of the current paper to construct 
similar advanced models for NGC 4494, but we do explore 
a flattened scenario using similar Jeans machinery to H4.2I 
as described in Appendix [C] We find that the inferred DM 
halo differs little from our spherical conclusion. It will be the 
object of a future paper to fully explore multiple modelling 
methods with NGC 3379, NGC 4494 and other galaxies, to 
get a handle on the robustness of the results to the methods. 
But so far, although very general models allow a wider range 
of solutions, we have seen no evidence that these models 
imply a systematic shift in the preferred mass solutions. 

Fully independent mass constraints from another obser- 
vatio nal tracer such as GCs or X-r ay emitting gas are desir- 
able. lO'Sullivan fc Ponmanl (|2004T) analysed XMM-Newton 
observations of NGC 4494, finding a weak component of very 
soft (fcaTx ~ 0.2 keV) X-ray emitting gas in a compact ar- 
rangement (within ~ R c ). They proposed that the X-ray 
underluminosity of this and similar ETGs could be due to 
either a weak or diffuse DM halo, or to a young field star 
population which has had little time to build up an X-ray 
halo through winds (though there are problems with this 
scenario) . More recently, it was suggested that the apparent 
emission from hot gas in such galaxies may in fact be domi- 
nated by cataclysmic variables and coronally active binaries 
i|Sazonov et al.lbood : iRevnivtsev et alj|2008h . 

Whether or not it is feasible to measure the halo mass 
of NGC 4494 using hot gas is thus a very open ques- 
tion. Assuming that there is any such gas which could be 
studied with deeper observations, it is more likely to be 
in an outflowing wind state than in hydrostatic equilib- 
rium, which would grea tly complicate the mass analysis 
(Pellegrin i fc CiottiluOOq ). Given these caveats, we note that 
iFukazawa et al.l (|2006T ) did derive a basic mass estimate for 
NGC 4494 from Chandra observations of apparent compact 
X-ray gas emission. They found an M/L of 6.2 ± 1.9 Tb,q 
inside 1 R E (converted to our distance, with minimal un- 
certainties estimated from the uncertainty in Tx), which 
is consistent with our finding of ~ 4.5 T_g,© at this ra- 
dius (see Fig. 1 11 ft - Not e though that a different analysis by 
iDiehl fc Statlerri|2007f) found no trace of X-ray gas in this 
system. 



5 IMPLICATIONS: ORDINARY ELLIPTICALS 
REVISITED 

Given the dynamical solutions found for NGC 4494 in 34.21 
we compare them to other galaxies, and consider the possible 
implications for DM and for galaxy formation. We start with 
a small sample of galaxies consisting of all four "ordinary" 
ellipticals discussed in R+03, whose PN dynamics origi- 
nally suggested unexpectedly low DM content: NGC 821, 
NGC 3379, NGC 4494 and NGC 4697 (with data from the 
latter presented in iMendez et afl feoOl). These galaxies have 
similar central velocity dispersions ctq ~ 200 km s~ x and 
near-L* luminosities (Mb ~ —20.3; see Table|6]), and are of 
the discy/cuspy ETG sub- type found by N+05 to have lower 
DM content on ~ 5R e scales than the boxy/cored ETGs. 
These galaxies (dubbed the "R+03 sample") are also of 
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Table 6. Mass results for fast-rotator early- type galaxies, based 
on extended kinematics data and higher-order LOSVD modelling. 
The absolute magnitudes Mg and central velocity dispersions 
Mg and ctq are adapted from HyperLeda, D+07, and this work. 



Name 


M B 


cro ( km s 1 ) 


log M« 


log M vir 


c vir 


Tvir 








(M ) 


(M ) 




(TTq.b) 


NGC 821 


-20.5 


210 


11.24 


14.46 


2.5 


1700 


NGC 3379 


-19.8 


207 


11.01 


12.17 


6.8 


110 


NGC 4494 


-20.5 


150 


10.99 


12.06 


7.6 


47 


NGC 4697 


-20.2 


174 


11.14 


12.80 


6.6 


300 



the "fast rotator" variety of ETG ( H3.2l and Emscllcm et al.l 
120071 ). which as a class shows hints of less central DM con- 
tent than the slow rotators (|Cappellari et al1l2006l ). 

Each galaxy has now been modelled using constraints 
on higher-order LOSVD moments to handle the mass- 
anisotropy degeneracy, and to infer the T* component by dy- 
namical means: NGC 821 using three- integral axisymmetric 
orbit models (stellar kinematics only; iForestell fc Gebhardtl 
120081 . but see also prel i minar y analysis of PN kinemat- 
ics in iNapolitano et alj 120071 ): NGC 3379 using spheri- 
cal orbit models (R+03) and spherical, axisymmetric and 
quasi-triaxial particle modelling (DL+08b); NGC 4494 us- 
ing kurtosis-based Jeans models; and NGC 4697 using ax- 
isymmetric particle modelling (DL+08a) . One caveat is that 
these results are tied to particular assumptions about the 
form of the DM density profile (NFW or LOG), which affects 
the decomposition of the total mass into stars and DM. E.g. 
dynamics on its own does not completely rule out a scenario 
(though unrealistic) where the stars have zero mass and the 
central DM profile is similar to the stellar luminosity profile. 

We examine the DM fraction for these galaxies, as well 
as for theoretical models, in H5.ll In H5.2I we consider the 
global DM halo parameters, and in H5.3I we place them in 
the wider context of galaxy formation. 

5.1 The dark matter fraction 

The first basic parameter we consider is the DM fraction 
within a given radius r: 



/dm (r) 



M DM (r) 



(17) 



Mtot(r) ' 

where we find for NGC 4494 from the best-fit NFW and 
LOG models that at 5 R c , /dm, 5 = 1 — T*/Ts ~ 0.2- 
0.5. For NGC 3379, the comparable calculation from the 
R+03 spherical model yields / DM , 5 ~ 0.2-0.4 (cf. D+07). 
The flattened models for NGC 3379 from DL+08b imply 
/dm,s ~ 0.2-0.6, while similar models for NGC 4697 give 
/dm, 5 ~ 0.3-0.6. These galaxies are therefore all consistent 
with a typical value of /dm ~ 0.4 at 5 R e ; this value is 
notably lower than one would expect from strong gravi- 
tational lensing studies whi ch typically find /dm ~ 0.3 at 
R e (e.g. iGavazzi et al.ll2007r i. i.e. at much smaller distances 
from the centre. However, lensing studies so far include very 
few of the fainter ordinary galaxies which would be equiva- 
lent to this fast rotator sample, altho ugh there are emerg- 
ing indications of lower /dm systems (|Ferreras et alj|2005l ; 



ICovone et ai1l2008l ). The constraints on NGC 821 (based on 
stellar kinematics only) do not extend to large enough radii 
for a direct comparison, and the models do suggest that this 
galaxy is different, with /dm ~ 0.5 at 1.5-2 R c already. 
However, our prelimina ry analysis based on PN kinematics 
ijNapolitano et al1l2007h seems to favor a lower /dm ~ 0.4 
at 5 R a , as found for the other three observed galaxies. 

On the theoretical side, N+05 found for a simpli- 
fied ACDM-based suite of models that one should expect 
/dm, 5 = 0.5-0.8, depending on various assumptions about 
the galaxy and DM masses. A similar framework from ML05 
predicted /dm, 5 ~ 0.5. Simulations of ETG formation in- 
cluding baryonic processes are now beginning to produce 
plausible results, with the schematic pair-merger simulations 
of D+05 yiel ding /dm. 5 ~ 0.5- 0.7. The full cosmo l ogica l 
simulations of iNaab et al.l (|2007h and lOfiorbe et all (|2007h 
suggest /dm, 5 ~ 0.4 and ~ 0.6-0.7, respectively (though it 
is not clear if these simulated galaxies are better analogues 
of the fast or slow rotators). This variety of predictions are 
all roughly consistent with a typical /dm, 5 ~ 0.5 (a quantity 
that is presumably not very sensitive to baryonic effects) — 
which appears similar to the empirically-obtained values. 

For additional insight, as in D+07 we show the DM 
fraction as a function of effective radius in Fig. 1151 with var- 
ious of the aforementioned empirical and theoretical results 
displayed. The plot is somewhat tricky to interpret because 
the stellar and DM mass profiles probably do not scale with 
luminosity i n a h omologous way. E.g. NGC 3379 and the 
INaab et al.l l|2007l ) simulated ETGs have similarly small val- 
ues for R e and hence should be compared; the larger galax- 
ies NGC 4494 and NGC 4697 shoul d presu mably be com- 
pared to the D+05 and Na ab et alj (|2007T l simulations. It 
thus becomes clear that simulations which include baryonic 
effects predict systematically more DM within 5 R e than 
is inferred observationally, and that this mismatch becomes 
even stronger at smaller radiF 3 ! 

While the dark matter fraction /dm, 5 is a directly- 
constrained observational quantity, its meaning may be 
muddled by the whims of the baryonic processes in in- 
dividual galaxies. For example, NGC 4494 appears to be 
more DM dominated in its central regions than NGC 3379 
(Fig. I15[) . but the global halo properties of the two galaxies 
turn out to be very similar — even though the PN velocity 
dispersion in NGC 3379 declines much more steeply than 
in NGC 4494 f ^O)) . This suggests that halo VD gradients 
are strongly affected by the individual galaxy scale- lengths, 
and careful dynamical modelling is necessary to interpret 
the implications for DM. 



5.2 Halo parameters: the R+03 sample 

For less ambiguous comparisons to theory, one can instead 
consider directly the properties of the DM haloes. A halo's 



13 Note that the uncertainty bounds on the observational results 
are not shown in the Figure for the sake of clarity. On a case- 
by-case basis, few of the galaxies may be formally inconsistent 
with the theoretical predictions, but the overall mean pattern is 
a systematic offset between observations and theory — while there 
is no indication of a systematic shift in the observational results 
(see also i|5,2l l. 
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Figure 15. Cumulative dark matter fraction as a func- 
tion of radius. The green solid line shows the best-fit /9(r) 
stars+NFW model while the short-long-dashcd blue line is the 
/3(r) stars+LOG model for NGC 4494: typical uncertainties are 
20-30% at 5R C (see Table and drop to ~ 10% at 1 R c . Light 
green long-dashed and cyan dot-short-dashed lines are the best- 
fit results of NGC 3379: DL+08b (model C90 - lower line - and 
D90 - higher line -) and R+03, respectively. Orange dot-long- 
dashed is the dark matter fraction of NGC 4697 (DL+08a model 
G with errorbar showing the DM fraction at 5 -Re for models 
D to J). The purple line with surrounding shaded region shows 
the mean results and scatter from the galaxy merger simu l ations 
of D+05. The red lines show simulatio ns from lNaab et all (2007, 
solid curves) and lOnorbe et al ] l|2007l . dashed curve). The black 
error-bar shows the theoretical prediction from N+05 based on 
ACDM cosmology and the full range of allowed star formation 
efficiency (egp =0.07-1; see Sec. 15.21 . 



virial properties are not directly constrained by PN kine- 
matics data which extend to only ~ 0.1r v i r , but we can 
at least hazard some guesses by assuming a priori a func- 
tional form for the mass profile to be extrapolated outwards 
(e.g. the NFW model). Dynamical modelling results based 
on NFW halo es were presented for NGC 82 1, NGC 3379 and 
NGC 4494 in iForestell fc Gebhardtl (|2008l ). R+03 and this 
paper, respectively. For NGC 4697, DL+08a used a LOG po- 
tential (model H), to which we have fitted an NFW model 
post hoc. The results are summarized in Table [6] where it 
should be noted that a cluster mass halo is not credible for 
the isolated galaxy NGC 821, and is a product of extrapo- 
lating from data extending to only 2 R c . However, the joint 
mass-concentration result for this galaxy is relevant, as dis- 
cussed below. 

Disregarding NGC 821, we see that the other three 
galaxies have virial M/L values of ~ 50-150 Tq^, and 
virial DM fractions of ~91%-95 %. Assuming a primordial 
baryon fraction of 0.17 l|Hinshaw et al.l [20081 ). this implies 
that the net star formation efficiency was ~ 30%-50% for 



these galaxies. These results are broadly consistent with the 
current consensus that star formation was most efficient in 
galaxies near the L* luminosi t y (e.g.. | Marinoni fc Hudson! 
120021 : iMandelbaum et alj|2006l ; Ivan den Bosch et al.ll2007l ). 

Rather than measuring the virial mass, the kinemat- 
ics data are best suited to measuring the amount of DM 
within ~ 5 R c , and to some extent the detailed profile of DM 
with radius in the same region. In the context of an NFW 
model, we can determine jointly the mass and concentra- 
tion parameters M v ir and c v j r . Given perfect observational 
constraints from small to large radius, we could uniquely 
determine these parameters as well as T* based on the de- 
tailed shape of the v c (r) profile. With real- world uncertain- 
ties in the observations and the dynamical modelling, there 
is in practice a strong degeneracy between these parameters, 
within certain limits (cf. Fig. [Til and ML05 Fig. 1). Roughly 
speaking, a larger value of M vlr will require a smaller value of 
Cvir (in order to decrease the fraction of DM inside 5 R e ) and 
a higher value of Y„ (compensating for the reduced amount 
of DM in the very central regions). 

We show the mass-concentration results (Table for 
the R+03 sample in Fig. [16] (filled triangles). Except for the 
case of NGC 4494, the uncertainties are generally not shown 
because at this juncture they are not uniformly and ro- 
bustly determined. As previously discussed, the uncertainty 
regions can be expected to follow a diagonal degeneracy 
track from low-mass and high-concentration, to high-mass 
and low concentration. This suggests that the NGC 821 re- 
sult at M v i r ~ 10 14 Mq and c v ir ~ 3 could be consistent with 
the other three galaxies at M v i r ~ 1O 12 M0 and c v i r ~ 6. 

Next we can compare these empirical results to theoret- 
ical expectations. In a collisionless ACDM universe, a mean 
relation is expected between mass and concentration such 
as the following: 



r(M v 
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M vir 
T0 n M 



(18) 



where the relation has a 1 a scatter of 0.14 dex, and is valid 
for z = 0, n m = 0.3, Qa = -7, h = 0.7, A vir = 101 and 
era = 0.9 (jBullock et al.l l200ll : N+05). Although more re- 
cent work has revised this relation with higher-resolution 
simulations and u pdated cosmologies (e.g. iNeto et al.ll2007l ; 
iDuffv et al.ll2008l ). the differences below M vir ~ 10 14 M Q 
are only at the level of ~ 10%. The prediction is thus that 
higher-mass haloes have lower concentrations, which conve- 
niently coincides roughly with the direction of the uncer- 
tainties in the observational results, and means that we can 
make useful comparisons to theory despite the modelling 
degeneracies. 

As shown in Fig. 1161 the R+03 sample suggests a mass- 
concentration relation that is parallel to the theoretical ex- 
pectation, but offset to lower concentration values by a fac- 
tor of ~2. This discrepancy was also manifested in ^5.11 as 
a difference in DM fractions. One could perhaps systemat- 
ically increase the observed concentrations into agreement 
with theory, by postulating that all of the galaxies system- 
atically have very radial halo anisotropy (e.g. /3 > 0.8 at 
~ 7 R e ; DL+08b) which may be technically allowed by the 
data. However, this anisotropy is not consistent with the 
ACDM-based simulations of D+05, nor perhaps with any 
reasonable theory of galaxy formation: for this ETG sam- 
ple, it appears difficult to simultaneously fit both the mass 
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Figure 16. Dark matter halo virial mass and concentration pa- 
rameters. Triangles and boxes mark fast-rotator and slow-rotator 
ETGs, respectively, "x " and "+ " symbols mark discy /cusped and 
boxy / cored galaxies respectively, classified following [Faber et al.l 
ll 19971 1 and N+05. The filled symbols mark detailed ETG dynami- 
cal results using PNe and GCs (light error bars, where available), 
with solid and dashed heavy error bars showing the statistical and 
systematic ^4. 311 uncertainties for NGC 4494, The open symbols 
show the dynamics- based ETG results from N+05, with error bars 
in the upper right corner showing the typical uncertainties. The 
dashed line shows the mean result for X-ray bright groups and 
clusters, the dot-dashed line is an inference for late-type galaxies, 
and the dotted line is the trend from weak lcnsing of all types of 
galaxies and groups. The solid curve with shaded region shows a 
mean relation expected from ACDM, with its 1 a scatter. 



and anisotropy profiles of ACDM galaxies. However, since 
the anisotropy is not yet well constrained by observations 
or by theory, we will continue to focus here on the better- 
understood issue of mass profiles. 

For comparing with models parameterized by the scale 
radius r s and density p s (e.g. Fig. I14p , we find that eqn. 1181 
is equivalent to the following relation: 



Ps 
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where the scatter in p s at fixed r s is a factor of 1.8. 

5.3 Halo parameters: wider context 

The mass-concentration comparison will obviously need 
much more scrutiny using a large observational sample 
along with rigorous dynamical modelling — which is the 
goal of our ongoing PN.S survey of ETGs. But based on 
the preliminary trend that we see, we may start putting 
the results in context, starting with some comparisons 
to previous results for ETGs. For the sake of unifor- 
mity, we initially include only studies of ETGs using ex- 



tended kinematics data fitted to basic NFW-based models 
as in this paper. These galaxies are NGC 1399, M87, and 
NGC4636 — in all cases modelled with a combination of stel- 
lar and GC dynamics, tho ugh not in all cases directly con- 
straining the anisotro p y dRomanowsky fc Kochanekl 1200 ll ; 
ISchuberth et all l200d : iRichtler et alJ bfJOSh . Intriguinglv. 
these galaxies' haloes [filled boxes in Fig. \16} all have much 
higher concentrations than the R+03 sample (by a factor of 
~3-4), and are slightly higher than the theoretical expecta- 
tion, though consistent within typical modelling uncertain- 
ties. 

We do not believe this difference is a systematic mod- 
elling effect from using PN vs GC dynamics, as our pre- 
liminary PN.S work al so suggests a populati on of high- 
concentration galaxies dNapolitano et aLll2007h. T he ETG 
stellar dynamics analysis of iGerhard et alT l 200 ll ) found a 



main sequence of high DM halo densities consistent with the 
three high-concentration galaxies, along with a few lower- 
density galaxies similar to t he four low-concentra tion galax- 
ies. A related analysis from lThomas et al.l (|2007l ) also found 
a large scatter in DM halo densities for cluster ETGs. 

In order to increase the sample size, we next plot the 
results from N+05, who used dynamical results on ETGs 
from the literature to infer NFW-based mass-concentration 
parameters. We classify the galaxies as fast- or slow-rotators 
based on C+07 where possible, and otherwise using litera- 
ture results as described in H3.2I While these mass inferences 
are not as individually reliable as direct NFW-based mod- 
elling, we do not detect any systematic difference between 
the two techniques when comparing overlapping cases. Inter- 
estingly, the main effect of this supplemental sample is to fill 
in the "missing" population of normal-concentration haloes 
(Fig. [TBJ . Examples of such galaxies include NGC 4406, 
NGC 5128, NGC 5846 and NGC 6703. 

When compared to the theoretical prediction, the over- 
all ETG dynamics-based results suggest a normal popula- 
tion of haloes, with a systematic bias for fast-rotators to re- 
side in low-concentration haloes, and slow-rotators in high- 
concentrations. Although it is outside the scope of this paper 
to determine if this difference is more strongly linked to ro- 
tation, luminosity, isophote shape, environment, etc., Fig. 1161 
suggests that isophote shape is at least as important as ro- 
tation. The concentration gap seen in the initial sample may 
be due to selection effects, e.g. the systems with the most 
GCs amenable to dynamic al study would be s ystematically 
the most massive systems (Spitl er et al.ll2008l ). But we also 
note that given no theoretical prejudice, the data themselves 
would suggest a shotgun pattern of masses and concentra- 
tions, with no hint of a systematic decrease of concentration 
with mass (if anything, the data suggest an increase). 

Since the initial three high-concentration systems are 
galaxies at the centres of massive X-ray-bright groups and 
clusters, we compare the resul ts from a large sam ple of low- 
redshift groups and clusters (D uffy et al.ll21IHKi. where we 
have transformed the halo parameters from A v i r = 200 to 
A v i r = 101, and from z = 0.1 to z = by the trend c v i r ~ 
(1 + z)~ l . As shown in Fig. 1161 this relation is similar to 
(~15% higher than) the high-concentration dynamical sam- 
ple. Other studies using gravitational lensing have also found 
relatively high concentrations for clusters ( Hcnnaw i et al.l 
| 2007 j; IComerford fc Nataralanll2007i ; iBroadhurst fc Barkanal 
booa broadhurst et al.ll2008l ). 
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We next broaden the picture to include late-type galax- 
ies. There is an emerging consensus that such galax- 
ies near the V luminosity (including the Milky Way) 
have NFW-based concentrations that are unexpectedly 



2008 



low jKassin et al.ll2006|; iDutton et al, Il2007l ; IMcGaugh et al 
20071 ; iGnedin et al.ll2007l ; IXue et al.ll2008l ; see lMaccio et al 



for a differing vi ew). We take th e mass - velocity re- 
lation from Fig. 10 of IMcGaugh et~ai1 (|2007h . and adapt 
it for Fig. 1161 where we see that it is similar to the low- 
concentration ETGs (~25% lower). 

The overall impression we gain from these combined 
detailed studies of individual systems is of an "S" -shaped 
relation, where the galaxy-scale haloes follow a trend of 
relatively low concentrations, and there is a transition at 
Mvir ~ 5 x 1O 12 M0 to group haloes with higher concentra- 
tions. This transition is comparable to the DM "dichotomy" 
between low- and high-stellar-mass galaxies found by N+05, 
occurring at M» ~ 2 x 10 n M Q . 

As a final piece to the puzzle, we consider a recent 
weak-gravitational-lensing study of both late and early-type 
galax ies, from galaxy to cluster masses (|Mandelbaum et al.l 
2008). The mean trend is shown in Fig. [16] (after conver- 
sions for Avir and z) and is seen to roughly coincide with 
the apparent trends for late-types and fast-rotator ETGs, 
with no indication of the high-concentration trend found by 
dynamics and by X-rays. 

Although there are potential sources of systematic er- 
ror with any of these mass-measurement techniques, and 
some of the sample sizes are small, the observational evi- 
dence so far agrees on a main population of relatively low- 
concentration galaxy-scale DM haloes. At higher masses, the 
observational inconsistency might be resolved if there were 
a strong selection effect in the X-ray and dynamical stud- 
ies for high-concentration systems. Alternatively, keeping in 
mind that weak lensing probes the outer regions of haloes 
while the other techniques probe the central regions, the 
inner haloes of groups and clusters might deviate strongly 
from the assumed NFW profiles. 

We next consider some theoretical scenarios that might 
explain these observational constraints. The impact of 
baryons on their host DM haloes could be a significant 
factor; e.g. the observed concentration transition occurs at 
roughly the same mass scale where a change in the physics o f 
gas cooling and heating is expected (|Cattaneo et all 12006). 
As discussed in H4.3I the standard model for baryonic mod- 
ification of a collisionless DM halo involves adiabatic con- 
traction, which draws additional DM into the central re- 
gions. The implication for inferring unmodified NFW con- 
centrations would be to decrease them, i.e. the observa- 
tional data points in Fig. [T5] would be shifted slightly down 
and to the left. This effect might bring the slow-rotator 
group-central ETGs into agreement with theory, but would 
make matters worse for the fast-rotator ETGs and the late- 
types. For those, a strong net halo expansion would be re- 
quire d, perhaps as a co n sequence of baryonic mas s expulsion 
(e.g. iMo fc Mad |2004 iMashchenko et all 120081 ). Alterna- 
tively, other effects such as dyna mical friction might flatte n 
the very central DM cusp (e.g. iRomano-Diaz et al.| [2008). 
which would not affect the overall concentration but could 
conceivably bias the observational inferences when fitting 
NFW mass models. However, the low concentrations found 



by weak lensing could not be ascribed to baryonic effects, so 
another explanation is needed. 

On a related note, it should be remembered that the 
NFW profile is not a prediction for a real halo even in the 
absence of baryons, as it represents a composite average of 
many haloes. Individual haloes would have departures from 
the NFW density, which might systematically connect to 
galaxy or group type in a way that creates the appearance 
of systematic concentration differences. 

As already mentioned, another possibility is that Fig. 1161 
reflects selection effects, i.e. the X-ray and GC studies are 
preferentially picking up systems that have high central DM 
concentrations (producin g unusually high X -ray luminosities 
and rich GC systems; e.g. lFedeli et al.ll2007l ). Because higher 
DM densities should arise from halo collapse at earlier times 
(when the background densities were higher), the implica- 
tion is that the high-concentration groups collapsed earlier. 
In order to not contradict the generic expectation of hierar- 
chical structure formation that higher-mass haloes collapse 
later on average (producing trends such as Eq. ll8p . the small 
observed subset of high-concentration systems would come 
from a high-cr tail of earlier-collapsing high-mass haloes. The 
different collapse times might also couple to baryonic effects, 
e.g. adiabatic contraction was stronger at earlier times when 
the gas content was higher. Such a scenario could be related 
to suggestions that the slow rotators as a class formed in 
rapid multiple mergers at higher redshifts, while fast rota- 
tors (and presumably late-types) formed at relatively low 
redshifts, experiencing a sm all number of isolated major 
mergers l|Burkert et al.ll2007l) . 

Thus, there may be a combination of effects in opera- 
tion: a mass dependency of baryonic, collapse, and merger 
processes in galaxy halo formation, coupled with an obser- 
vational bias. Whether such a combination could quantita- 
tively explain the observed concentration trends is a ques- 
tion far beyond the scope of this paper, but we do highlight 
some potential problems. 

One issue is the sheer magnitude of the difference be- 
tween the high- and low-concentration populations. The col- 
lapse redshift z c of a DM halo scales roughly with its cen - 
tral density p s as (1 + z c ) oc pl^ 3 (e.g.. iBullock et aill200ll ). 
The difference between the halo densities on the upper- and 
lower-concentration sequences is a factor of ~ 20, which im- 
plies markedly different formation redshifts at the same halo 
mass. E.g. if a typical 1O 13 M0 halo assembled at z c ~ 0.9 
(7 Gyr look-back time), then its high-concentration counter- 
part formed at z c ~ 4 (12 Gyr halo age). These haloes would 
comprise a very small fraction of th e total population of the 
same mass (e.g. lHarker e t al. 2006), whose presence should 
be difficult to miss by the dynamical and X-ray observations. 

The other problem is that the observed concentra- 
tions for the main halo population are too low. Unless 
there is a vast population of higher-concentration "dark" 
galaxy-scale haloes that are still missed by observational 
surveys, the remaining avenue is alteration of the cosmo- 
logical parameters. The obvious candidate, as, would have 
to be lowered substantially to produc e the low-concentratio n 
trend, perhaps to implausible values (jMcGaugh et al.ll2007l ). 
An alternative is to fine- tune the mass dependence of the 
DM power spectrum (e. g. Alam et al.|l2002 ■ IMcGaugh et al.l 



JJM power spectrum (e.g. Alam et al. Avva 
120031 ; iFedeli et al.ll2008l : iMaccio et alj f2008) 



The foregoing discussion is a speculative attempt to ex- 
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plain the entire observational picture in the context of the 
ACDM paradigm. However, there is still little direct evi- 
dence for ACDM halo profiles among all the observational 
studies of galactic and super-galactic mass distributions; e.g. 
in NGC 4494, the NFW and LOG potentials fit the data 
equally well. The possibility remains that the universe is 
dominated by an alternative form of DM or of gravity, which 
might more naturally explain the observations summarized 

in Fig. ma 



6 CONCLUSIONS 

We have presented observations of 255 PNe with line-of-sight 
velocities in NGC 4494, out to ~ 7 R e - The basic spatial and 
kinematical properties of these PNe (including rotation, dis- 
persion and kurtosis) agree well with the field stars, and do 
not support a distinct origin for the PN progenitors. The 
mean rotation is low but the system may become rotation- 
ally dominated in its outer parts, and we also see evidence 
for a kinematically-decoupled halo which may be a merger 
signature akin to kinematically-decoupled cores. The pro- 
jected velocity dispersion profile declines with the radius, 
although not as steeply as in NGC 3379. 

We have constructed spherical dynamical models of 
the system, including a pseudo-inversion Jeans mass model 
as well as a more standard Jeans analysis using multi- 
component mass models that include NFW or logarithmic 
DM haloes, with fourth-order moments being used to con- 
strain the orbital anisotropy. Both approaches give similar 
results for the mass profile within the radial range con- 
strained by the data. Some DM is required by the data; 
our best-fit solution has a radially anisotropic stellar halo, 
a plausible stellar M/L, and a DM halo with a fairly low 
concentration. 

We review the halo parameters determined observation- 
ally for other ordinary early-type galaxies such as NGC 4494, 
as well as for galaxies and groups from the literature. There 
are preliminary indications that most haloes follow a sim- 
ilar mass-concentration sequence which implies rather low 
DM densities, accompanied by a population of bright slow- 
rotator ETGs in group-scale haloes with much higher DM 
densities. We discuss some possible theoretical scenarios for 
tying these observations together. Clarifying the situation 
will require a broad observational survey of galaxies of dif- 
ferent types and environments, combined with rigorous dy- 
namical models and improved simulations of galaxy for- 
mation. It is the f ocus of our ongoing PN.S survey (e.g. 
ICoccato et"al]|2008l ) to fill in the critical observational gap of 
the halo properties of ordinary early-type galaxies. The pri- 
mary programme includes 12 ellipticals representing a range 
of sub-types and environments, with an extended survey of 
40 early-types in progress, limited only by volume and mag- 
nitude. 
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APPENDIX A: PHOTOMETRY OF NGC 4494 

We have obtained a composite photometric profile of 
NGC 4494 by combining new Megacam observations, de- 
tailed in H3.ll with literature data w hich include HST- based 
observations in the V and / bands (lLauer et al. 20051). and 
groun d-based CCD observations in BVI ( Goudfrooii et ail 
11994 hereafter G+94). The HST data are used inside 4.3", 
then G+94 to 32", and Megacam outside 32". The a 4 
isophote shape parameter l|Bender et al.|[l988l ) is not avail- 
able for the HST data. A small color gradient is found be- 
tween the G+94 and the Megacam data at the outer data 
points, probably due to the background subtraction. 

We derive some global photometric parameters for the 
galaxy by averaging the radially-binned values outside 5" 
(the central disc region), and weighted by bin luminosities. 
Thus we find e = 0.162 ± 0.001, al = 0.16 ± 0.07 and PA = 
-0.9° ±0.3°. 



APPENDIX B: JEANS EQUATIONS 

Here we present the second- and fourth-moment Jeans equa- 
tions and projections to velocity dispersion ai os and kurtosi s 
Kios, as needed in £14.2.31 and closely following lLokasI ([2002). 
It is the first time this particular method has been applied 
to an elliptical galaxy. 

We assume that the system is spherically symmetric and 
that there are no net streaming motions (e.g. no rotation) 
so that the odd velocity moments vanish and the different 
component of the dispersion velocity tensor are of = v%, 
(Tg — Vg and 00 = Vj,, with <rf = cr^ by symmetry. Under 
these assumption th e Jeans equations reduc e to the simple 
radial equation (e.g. iBinnev fc Mamon|[l982h : 



d f . 2 , , 2/3 . 2 .d* 



(Bl) 



where j*(r) is the 3D number density distribution of the 
tracer population (e.g. the stars or PNe) and <&(r) is the 
gravitational potential GM(r)/r 2 with G the gravitational 
constant. In equation (jBll) . the anisotropy parameter, (3 = 
1 — <jg /of , can be either constant or variable with the radius 
r. In the first case the equation is simply written as 
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Table Al. Combined V-band surface photometry of NGC 4494 
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second case the general solution is 
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00. In the 

(B3) 



under the same conditions. 

The line-of-sight velocity dispersion is obtained from 
the 3 D velocity dispersion by integrating along the line of 
sight (|Binnev fc Mamon|[l98a) 



I(R) 



R 2 



3* °l r 
Vr 2 - R 2 



dr 



(B4) 



where I(R) is the surface density of the tracer and R is 
the projected radius. For various simple functional forms for 
/3(r) (e.g. P constant), the calculation of ai os can be reduced 
to a one-dimensional numerical integration: 



o-ioa(-R) = 



2G 
1(B) 



Kj*Mdr 



(B5) 



where K(r) is a special function depending on /3(r) (see 
Appendix A of ML05). 

To address the well-known mass-anisotropy degeneracy, 
we consider higher-order moments of the velocity distribu- 
tion. For the fourth-order moments, the three distinct com- 
ponents Vr, v\ — vt and v 2 v 2 = v?v 2 are related by two 



higher order Jeans equations (|Merrifield fc KentHl99(j i 

In order to solve these equations we need additional in- 
formation about the distribution function. We therefore re- 
strict ourselves here to functions which can be constructed 
from the energy-dependent distribution function by multi- 
plying it by a function of angular momentum f(E,L) = 
fo(E) L~ 213 with P = const . This is an ansatz wide! 



iy 

used (Henonll973l. Dcjonghc 1986, Wil kinson fc Evan J 19991 . 
I An fc Evandl2006r ). which has the advantage of being easy 
to integrate even though it does not allow to generalize to 
the case of /3 — f3(r) for the 4 th -order moment. 
The two Jeans equations then reduce to one: 



dr r 
whose solution is 
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ijlokasl 120021 ). By projection, we obtain the line-of-sight 
fourth moment 
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which we calculate numerically. 

The fourth projected moment is used together with af os 
in order to obtain the reduced projected kurtosis: 



Klos(-R) = 



< S (R) 



(B9) 



which takes the value of for a Gaussian distribution. The 
advantage of using the kurtosis is that it can be easily de- 
rived from a discrete radial velocity distribution such as our 
PN data set (see Section \3. 41) . 
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One useful corollary of this formalism is that if <7i os is 
observationally constant with the radius R, and one assumes 
that the anisotropy is constant, then the internal kurtosis 
n(r) = 0, and one can calculate a one-to-one relation be- 
tween Kios and p. Therefore one can estimate the anisotropy 
as a direct deprojection of the observations, with no dy- 
namical modelling necessary. This approach is of course not 
applicable to NGC 4494 with its declining VD profile. For 
reference, the formulae are the following: 



Klos(-R) = 3 ( 
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APPENDIX C: FLATTENED MODEL 

Here we present an alternative model needed for i]4.3l 
wherein the galaxy is assumed to be axisymmetric and in- 
trinsically flattened but apparently round because it is seen 
near face-on. To simplify the analysis, we study the dynam- 
ics along the galaxy's observed minor axis only, where the 
motions in the galaxy's equatorial angular direction (j> do 
not contribute to the observed line-of-sight velocity disper- 
sion <Tl os . 

Our starting point for the Jeans analysis is an approx- 
imate relation between the vertical velocity dispersion and 
potential in a flattened galaxy: 



o z (R,Z) = - 
J* 



*dz dZ > 



(CI) 



where 7? is the equatorial radius a nd Z is the vertical height 
in the galaxy's reference frame fcf. lBinnev fc Tremainef"l987l . 
eqn. 4-36). Since we are observing along the minor axis, 
we will not be concerned with the azimuthal dispersion o$, 
while for the radial dispersion or we adopt an anisotropy 
parametrisation: Prz = 1 — o~ z/&r- 

The local contribution to the dispersion along the line- 
of-sight 2 is: 

°" 2 = °\ shi 2 i sin 2 (f> + cr 2 , sin 2 i cos 2 (f> + ct| cos 2 i, (C2) 

where i is the galaxy's inclination relativ e to z (see Fig. IC1I 
and c.f. iRomanowskv fc Kocha nck 1997 eqn. 17). Along the 
minor axis, <j> = 90° and we have: 



2 2 • 2 • 2 2 ■ 2 

°~z = °~R sm Z + C" z COS I = O z 



1- Pb 



+ cos i 



(C3) 



The observed dispersion in the sky plane (x,y) is then cal- 
culated by: 



o-\ as {x,y) 



j t o z dz. 



(C4) 



With a luminosity model j* (R, Z) and some assumptions for 
i and /3rz, we now have sufflcent equations to link $ to o"i os . 



.a, 



/Of- 



observer 




Figure CI. Schematic geometry of the flattened galaxy model, 
viewed from a line-of-sight in the plane of the sky. Shown are the 
inclination i, the intrinsic major and minor axis lengths a and 
6, corresponding to radial and vertical velocity dispersions ur 
and oz- The galaxy is axisymmetric in the equatorial plane (i.e. 
appears circular when i = 0°). 



We derive the intrinsic ellipticity £30 from the observed 
ellipticity £2D and an assumed inclination i using the follow- 
ing equation: 

(1 — 6 2 d) 2 = (1 — £3d) 2 sin 2 i + cos 2 i (C5) 

ijBinnev fc Merrifieldl[l998l . eqn 4.29). For a razor-thin disc 
(e3D = 1), the minimum inclination allowed by the ob- 
served flattening of NGC 4494 (e 2D ~ 0.14 for R Z 0.1 R c ) 
is i — 31°. More realistic ally, the most flattened E TGs 
observed have e 2D ~ 0.75 (|Binnev fc Merrifieldl [l998h . so 
as an extreme alternative to the spherical case, we adopt 
e 3 D = 0.75 (b/a = 0.25) and i = 32°. For equation IC3l we 
approximate this to i = 30° and find: 



2 

o- z 



4(1 -Prz) 



+ 



(C6) 



We next assume the vertical anisotropy Prz to be con- 
stant throughout the galaxy, and therefore make use of the 
SAURON dynamical results for galaxies' central regions to 
make an a priori estimate for Prz- From C+07 Fig. 7 we 
find that Prz ~ 0.5 is typical for galaxies with the flatten- 
ing of NGC 4494 (e 3D = 0.75). From equation IC6l this im- 
plies o~ z — 1.1 x oz, i-e. the line-of-sight dispersion is closely 
related to the intrinsic vertical dispersion. The anisotropic 
flattening of the galaxy means oz is relatively low, which 
is exactly the scenario we are attempting to explore for 
explaining the low observed values of <ti os . Note also that 
higher values for Prz will for the same data imply lower 
intrinsic oz and thus lower mass since the contribution to 
Oz from or will increase. 

As a consistency check on our model assumptions, 
we consider the rotation dominance parameter within R e , 
(V/o) c . For the SAURON galaxies, this parameter corre- 
lates well with the ellipticity when both parameters are de- 
projected (C+07 Fig. 3). For an NGC 4494 observed value 
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of (V/a) e Ri 0.26 (fj372j , we search for values of i and (5rz 
that bring the galaxy into agreement with the correlation 
(cf. C+07 eq. 12). We find i ~ 40° and /3 RZ ~ 0.2 with in- 
trinsic £3d ~ 0.4 and (V/a) c ~ 0.4, suggesting that the true 
parameters of NGC 4494 could well be intermediate to the 
extreme spherical and flattened cases that we are consider- 
ing. Note also that advanced axisymmetric modelling of the 
comparable galaxies NGC 3379 and NGC 4697 found that 
(3rz increases fr om ~ 0.3 in the central regions to ~ 0.5 in 
the outer parts (|De Lorenzi et al.ll2008al lbh. 

Next we follow the same steps as in the spherical case 
f £|4.2.3[) to explore a series of galaxy+DM mass models and 
c ompare them to t he data. For j*(z) in Ean. [C4l we adopt 
a iHernquistl (|l990h profile for simplicity with the approxi- 
mation tha t the Hernqu ist radius is the usual th = 0.55 R c 
and r — y/ R 2 + (z/q) 2 , with q — 0.25, in order to account 
for the flattening along the line-of-sight. 

Recall that we are now comparing to « rms data along 
the minor axis only, with no attempt to fit the kurtosis. The 
best-fit dark matter halo has T* = 3.9X0,3, logMoM = 
11.91±0.20 (the dark mass), a concentration parameter c<j = 
2 ± 1, and a total virial ratio / V i r = 8 ± 5. The latter is lower 
than the lowest value obtained for the spherical case. As 
a comparison, we found T(i? v i r ) = 35 and Tb5 — 5 ± 1, 
much lower than the values found in Table [S] We therefore 
find that a face-on flattened configuration is not sufficient 
to produce a very massive DM halo solution for NGC 4494. 
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